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ABSTRACT 


In  salvo  firing,  the  smallest  ballistic  dispersion  is 
not  always  most  desirable.  Deliberate  increase  of  the 
ballistic  dispersion  can  improve  the  probability  of  des¬ 
troying  the  target.  Our  concern  in  this  dissertation  is 
the  optimization  of  such  "artificial"  dispersion  in  two- 
dimensional  salvo  models.  In  some  cases  no  closed  form 
solution  is  available,  but  we  are  able  to  offer  efficient 
methods  for  the  computation  or  approximation  of  the  salvo 
kill  probability.  In  other  cases  we  are  able  to  derive 
approximate  formulae  for  the  optimal  ballistic  dispersion 
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I .  INTRODUCTION 


Let  us  begin  with  a  classical  example  presented  in  Chap¬ 
ter  6  of  the  textbook  by  Morse  and  Kimball  [Ref.  1] .  Suppose 
an  airplane  carries  two  bombs  to  attack  a  railroad  track. 

The  plane  flies  along  a  course  perpendicular  to  the  track, 
and  drops  the  bombs.  There  are  three  ways  to  drop  two 
bombs:  (1)  together  in  salvo,  (2)  spaced  a  suitable  distance 

apart  aiming  the  midpoint  of  the  pattern  at  the  center  of 
the  track,  and  (3)  dropping  each  on  a  separate  run  over  the 
target . 

We  intend  to  illustrate  the  method  of  calculating  the 
probability  of  track  destruction  in  the  three  cases,  the 
purpose  being  to  determine  which  method  of  attack  is  the 
best. 

Let  us  consider  Case  (1) ,  salvo  bombing.  The  two  bombs 
leaving  the  plane  simultaneously  will  hit  the  ground  at  some 
distance  apart;  the  impact  point  being  random.  For  the  time 
being  the  random  vector  of  the  impact  position  from  the 
center  of  impact  is  termed  the  ballistic  error.  The  aiming 
is  also  not  free  from  error.  The  plane  aims  at  a  point  on 
the  center  of  the  track,  but  the  center  of  impact  will 
deviate  from  the  aim  point.  This  deviation,  or  aiming  error, 
is  also  random.  The  impact  point  of  each  bomb  is,  there¬ 
fore,  composed  of  two  errors;  first  the  aiming  error  which 
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is  common  to  both  bombs,  and  second  the  ballistic  error 
which  varies  from  one  bomb  to  another. 

If  the  variance  of  the  ballistic  error  is  small  while 
that  of  the  aiming  error  is  large,  the  two  bombs  will  land 
almost  at  the  same  point.  Therefore  it  is  expected  that 
either  both  will  hit  or  neither  will  hit  the  track.  To  avoid 
the  latter  situation  and  to  improve  the  probability  of  des¬ 
truction,  it  is  better  to  spread  the  landing  pattern  of  the 
bombs.  One  way  is  the  method  of  pattern  bombing  described 
as  Case  (2) ,  another  is  salvo  bombing  employing  bombs 
with  a  suitable  value  of  variance  of  the  ballistic  error. 

Our  main  concern  in  this  dissertation  is  maximization  of  the 
probability  of  salvo  destruction  by  choosing  a  suitable 
value  of  the  variance  of  the  ballistic  error.  In  the  follow¬ 
ing  we  call  it  the  optimal  ballistic  dispersion. 

The  problem  described  above  is  not  a  new  one.  In  an 
interesting  talk  before  a  conference  at  the  Ballistic  Re¬ 
search  Laboratories  held  in  March  1955  iRef .  2J ,  Merritt 
and  King  stated  that  an  approximate  formula  for  the  optimal 
ballistic  dispersion  was  derived  by  two  Englishmen  as  early 

as  in  1936.  Following  the  work  in  WW  II,  a  number  of  arti- 

\ 

cles  were  published  on  both  the  salvo  and  pattern  firing 
models,  but  since  the  early  1970's,  it  seems  that  these 
models  have  attracted  less  attention.  Further  detail  can  be 
found  in  two  excellent  review  papers  by  Eckler  [Ref.  3]  and 
Eckler  and  Burr  [Ref.  4] .  However,  little  is  known  as  to  the 


optimization  of  the  ballistic  dispersion  in  salvo  models. 

It  is  the  calculation  or  approximation  of  the  optimal 
ballistic  dispersion  that  is  the  subject  of  this  dissertation. 

A.  SALVO  FIRING 

In  the  following  chapters  we  deal  with  only  two-dimensional 
salvo  models.  The  reason  we  adopt  the  two-dimensional  model 
is  that  it  is  the  most  frequent  case  and  plays  an  important 
role  in  real  world  applications.  One-  or  three-dimensional 
theory  can  also  be  developed  to  parallel  our  investigation. 

Suppose  that  there  is  a  target  in  a  two-dimensional  space, 
and  a  salvo  of  n  weapons  is  delivered  against  the  target. 
Delivery  error  relative  to  the  target  is  assumed  to  be  com¬ 
posed  of  two  parts,  the  aiming  error  and  the  ballistic  error. 

First  we  aim  at  the  target.  With  this  aiming,  the  center 
of  impact  point  of  n  weapons  is  determined.  Let  us  adopt  a 
Cartesian  coordinate  system  (x,y)  such  that  the  origin  coin¬ 
cides  with  the  center  of  impact.  The  aiming  cannot  be  per¬ 
fect.  Let  the  components  of  the  aiming  error  be  -U  and  -V. 

Then  the  position  of  the  target  with  respect  to  our  coordinate 
system  is  given  by  (U,V) .  It  is  assumed  that  the  joint 
probability  density  of  U  and  V  exists  and  is  given  by  f1(u,v) . 

Now,  a  salvo  of  n  weapons  is  fired  against  the  target 
after  (U,V)  takes  a  value  (u,v)  which  is  unknown  to  us  and 
only  is  predictable  in  a  probabilistic  sense.  The  components 
of  the  impact  point  of  the  ith  weapon  are  denoted  as  X^  and 
Y..  It  is  assumed  that  (X.,Y.)  are  independent  and  identically 


distributed  random  variables  with  joint  density  f2(x,y) . 

In  Fig.  1,  a  typical  geometry  is  shown.  The  target  is  at 
(u,v) ,  and  the  impact  points  of  four  weapons  are  scattered 
around  the  origin. 

As  to  the  target  we  adopt  the  so-called  point  target 
concepts  It  is  assumed  that  the  target  is  either  completely 
destroyed  or  else  undamaged  by  each  weapon.  We  neglect  any 
possible  partial  damage  and  its  cumulative  effect.  The 
probability  that  a  weapon  landing  at  (x,y)  destroys  the  tar¬ 
get  at  (u,v)  is  a  function  of  u-x  and  v-y,  denoted  as 
D(u-x,v-y)  and  is  called  the  damage  function. 

The  probability  of  destroying  the  target  is  calculated 
in  two  steps.  First,  suppose  that  the  random  variables 
U  and  V  take  some  values  u  and  v.  The  conditional  proba¬ 
bility  that  the  ith  weapon  destroys  the  target  given  U  -  v 
and  V  =  v  is 

oo  00 

P ( u, v)  =  /  /  D( u-x, v-y) f 5 (x,y) dxdy  .  Cl.l) 

— 00  —00 

The  impact  points  and  are  all  assumed  to  be  independent 
so  the  conditional  survival  probability  of  the  target  given 
U  =  u  and  V  =  v  is  {1  -  p(u,v)}n.  Therefore,  we  obtain  the 
probability  that  the  target  is  destroyed  by  the  salvo  of  n 
weapons— salvo  kill  probability — by  averaging  the  conditional 
survival  probability  with  respect  to  the  distribution  of  u 


and  v,  and  then  subtracting  it  from  unity: 


CO  00 

p  “  1  -  /  /  {1-  p(u,v) }nf1 (u,v) dudv  .  (1.2) 

.CO  —00 

B.  UPPER  BOUND  ON  THE  SALVO  KILL  PROBABILITY 

In  the  example  mentioned  earlier,  Morse  and  Kimball  ob¬ 
served  that  the  probability  of  destroying  the  track  by  drop¬ 
ping  two  bombs  in  salvo  is  always  smaller  than  the  destruc¬ 
tion  probability  of  Case  (3) ,  namely  the  case  in  which 
the  two  bombs  are  dropped  independently  on  separate  passes. 
Similar  observations  are  also  pointed  out  in  the  review  paper 
by  Eckler  and  Burr.  The  observation  is  valuable,  but  there 
has  been  no  proof  given  to  this  fundamental  property  dis¬ 
covered  numerically.  Indeed,  we  have  the  following  proposi¬ 
tion  without  any  specific  assumptions  on  f1(u,v),  f2(x,y) 
and  D(x,y) . 

Proposition  1.1 

Let  the  salvo  kill  probability  with  n  weapons  be  denoted 


n 


=  1  -  /  /  {1  -  p(u,v) }nf1 (u,v) dudv 


Then, 


(1.3) 


pn  <  1  -  (1-p,) 


(1.4) 


Namely,  salvo  Pr  is  always  smaller  than  the  kill  probability 
of  n  independent  tries. 


Proof 

The  proof  is  immediate  from  a  theorem  in  the  book  by 
Hardy  and  others  [Ref.  5] .  Define  a  mean  M  of  a  non-negative 
function  g(x)  for  a  positive  r  as. 


Mr(g)  =  {  /  grf  dx}1/r  , 


where  f  =  f  (x)  is  a  weighting  function,  positive  everywhere 


and, 


/  f  dx  ■  1  . 


Theorem  192  from  Ref.  5  tells  us  that  if  0  <  r  <  s  and 

M  (g)  is  finite,  then 
s 


M  <  M  , 
r  s 


unless  g(x)  is  a  constant. 

The  weighting  function  f  corresponds  to  f^(u,vl  in  our 
proposition,  and  g  corresponds  to  our  l-p(u,v).  Therefore, 
for  n  >  1,  we  have 
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Mj^d-p)  <  Mn(l-p) 


or 


{M1(l-p)}n  <  {f^d-p)}11  , 


which  means 


d-PjT1  < 


1  -  P_ 


1  n 


Here  is  an  outline  of  the  succeeding  chapters:  In  Chap¬ 
ter  II,  we  investigate  the  simplest  salvo  model  which  is 
characterized  by  circular  normal  errors  and  circular  Gaussian 
damage  function.  After  that,  an  approximate  formula  for  the 
optimal  ballistic  dispersion  is  introduced  in  Chapter  III; 
we  regard  the  material  in  Chapter  III  as  the  "core"  of  this 
dissertation.  In  Chapter  IV,  the  approximate  formula  is 
applied  to  the  salvo  model  of  Chapter  II,  and  its  accuracy 
is  studied.  In  Chapter  V,  an  approximate  method  of  calcu¬ 
lating  the  salvo  kill  probability  of  the  salvo  model  with 
the  general  damage  function  is  presented;  this  case  is  more 
difficult  than  the  one  where  the  damage  function  is  circular 
normal,  so  we  resort  to  further  approximations.  In  Chapter 
VI,  we  investigate  the  salvo  model  with  the  so-called  cookie 
cutter  damage  function,  relying  heavily  on  the  results  of 
Chapters  III  and  V.  Finally,  salvo  model  without  circular 
symmetry  is  discussed  in  Chapter  VII.  Our  concern  throughout 
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is  the  optimization  of  artificial  dispersion.  In  some  cases 
no  closed  form  solution  is  available,  so  we  offer  efficient 
methods  for  the  computation  or  approximation  of  salvo  kill 
probability  with  a  given  dispersion.  In  other  cases  we  are 
able  to  offer  closed  form  solutions. 


II.  SIMPLE  SALVO  MODEL 


In  this  chapter  we  deal  with  the  simplest  salvo  model 
with  circular  normal  errors  and  a  circularly  symmetric 
Gaussian  damage  function. 

Let  us  adopt  a  Cartesian  coordinate  system  (x,y)  the 
origin  of  which  coincides  with  the  center  of  impact  of  the 
n  weapons.  The  target  position  with  respect  to  this  coor¬ 
dinate  system  is  denoted  as  U  and  V  in  the  x  and  y  direction, 
respectively.  (U,V)  is  a  random  bias  common  to  all  the  n 
weapons.  In  this  chapter  (U,V)  is  assumed  as  a  circular  normal 
variate;  namely,  we  assume  that  U  and  V  are  independent  and 

identically  distributed  normal  random  variables  with  mean  0 
2 

and  variance  The  joint  density  of  U  and  V  is,  therefore 


f1(u,v) 


u2+v2 


2  a\ 


(2.1) 


The  impact  point  of  the  ith  weapon  is  denoted  as  (X^,Y^) . 
It  is  assumed  that  (X^,Y^) ,  i  =  l,2,...,n  are  independent  and 
identically  distributed  circular  normal  random  variables 
with  joint  density  function 


f2  (x,y) 


2.  2 


9- 


2a 


2ira! 


(2.2) 
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As  to  the  damage  function,  we  assume  that  the  so-called 


Gaussian  damage  function  with  circular  symmetry:  If  the 
target  is  at  (u,v)  and  a  weapon  impacts  at  (x,y) ,  then  the 
destruction  probability  is  given  by 

(u-x) 2+ (v-y) 2 
- 1  2“ 

D(  u-x,  v-y)  =  e  <£ct  .  (2.3) 

In  the  following,  let  us  call  this  model  the  Simple  Salvo 
Model. 

According  to  Grubbs  [Ref.  6],  the  Simple  Salvo  Model  was 
first  studied  as  early  as  in  1953  by  H.  K.  Weiss  in  BRL  Re¬ 
port  No.  879,  "Methods  for  Computing  the  Effectiveness  of 
Area  Weapons."  The  formula  for  the  salvo  kill  probability 
(2.7)  given  in  the  next  section  is  attributed  to  him.  Later 
Breaux  [Ref.  7]  found  that  this  formula  was  not  suitable  for 
computation  when  n  is  large,  and  gave  another  method  of  com¬ 
putation.  However,  it  seems  that  there  is  no  published  re¬ 
search  investigating  this  model  in  full  depth;  the  studies 
usually  come  to  an  end  when  expressions  for  the  kill  proba¬ 
bility  are  derived.  In  this  chapter,  we  investigate  first 
the  property  of  the  kill  probability  as  a  function  of  param¬ 
eters  involved  in  the  model.  Methods  for  computing  the  kill 
probability  and  bounds  for  the  kill  probability  are  also 
dealt  with  in  the  later  sections. 


A.  SALVO  KILL  PROBABILITY  AS  A  FUNCTION  OF  n,  c^,  o2  ,  AND  a 

The  salvo  kill  probability  is  derived  in  two  steps  as 
is  stated  in  Chapter  I.  First  the  conditional  kill  proba¬ 
bility  by  the  ith  weapon  is  calculated  given  that  the  random 
bias  (U,V)  takes  a  value  (u,v) . 


p(u,v)  *  /  /  D(u-x, v-y) f- (x,y) dxdy  . 


(1.1) 


— oo  —CO 


The  salvo  kill  probability  is  then  obtained  by  averaging 
1  - (l-p(u,v)  }n  over  the  distribution  of  (U,V)  ,  as 


p  =  1  -  /  /  (l-p(u,v) }nf1 (u,v) dudv  . 


(1.2) 


— OO  — oo 


Using  the  assumptions  (2.2)  and  (2.3),  we  easily  calcu¬ 
late  the  conditional  kill  probability  (1.1) 


2  2  2  2 
(u-x)  +(v-y)  x  +y 


and  get 


p(u,v) 


2  2 
u  +v 


2  2(a  +o:) 

a  2 


a  +  a. 


Substituting  (2.1)  and  (2.4)  into  (1.2),  the  salvo  kill 


probability  is 


-  1  - 


2tt  a! 


/  /  Cl  - 


a 

~T~~\ 

a  +a. 


2 ,  2 
u  +v 

- 2 — 

2  (a N-oJ) 


n 

}  e 


u2+v2 

2a? 


dudv 


Conversion  from  Cartesian  coordinates  into  polar  coordinates 
gives 


P 


a2 

~2  2 
a  +a2 


2  (a 


+02, 


n 


From  this  we  get 


=  1  -  /  (1  -  -  e"t/p)ne_tdt 


(2.5) 


where 


A  *  a 2/aJ,  p  =  (ct2+c^)/a^  .  (2.6) 

There  are  4  parameters  in  our  model,  but  the  salvo  kill 
probability  is  determined  by  three,  n,  o^/ot,  and  o2/ ra»  or 
n.  A,  and  p  as  is  seen  in  (2.5). 

If  the  integrand  in  (2.5)  is  expanded  in  a  binomial 
series  and  integrated  term  by  term,  we  obtain  the  Weiss 
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formula  for  the  salvo  kill  probability. 

P  -  l  (-D  (|) j  •  (2.7) 

j  =1  j  p  J  p 

P  is  obviously  an  increasing  function  of  n,  and  is  de- 
2  2 

creasing  in  o^/ot  .  The  latter  is  easily  verified  by  in- 

2  2 

specting  (2.5),  where  X/p  is  constant  for  fixed  a2/oi  *  Tlle 
2  2 

a2/a  -dependence  is,  however,  not  obvious.  In  (2.5)  ,  X  is 

2  2 

constant  for  a  fixed  value  of  c^/a  ,  and  there  are  two  factors 
which  involve  p  in  the  integrand,  with  p  tending  to  increase 
P  in  one  and  to  decrease  P  in  the  other. 

2  2 

Table  1  and  Table  2  are  given  here  to  show  the  a^/a  -  and 
2  2 

a2/ a  -dependence  of  the  salvo  kill  probability,  where  n  is 

kept  constant,  n  «  2  for  Table  1  and  n  =  16  for  Table  2.  For 

2  2  2  2 
fixed  n  and  a2/a  ,  P  is  decreasing  in  o^/a  as  mentioned 

2  2 

above.  However,  the  a2/a  -dependence  is  a  little  different. 

2  2 

It  is  observed  that,  generally,  P  is  decreasing  in  a2/a  for 

2  2  2  2 

small  values  of  a^/a  ,  but  for  a  larger  o^/a  ,  P  increases 
2  2 

when  o^/a  is  increased,  reaches  a  maximum  value,  and  then 

decreases.  It  is  also  noted  that  from  the  viewpoint  of  real 

2  2 

world  applications,  the  case  of  large  a^/a  with  large  n  is 

especially  interesting;  e.g.,  as  seen  in  Table  2,  the  kill 

probability  for  n  =  16  and  a£/a  =  8.0  is  improved  to  0.5881 
2  2 

when  a2/a  is  increased  to  a  suitable  value.  The  maximum  P 

2  2 

is  more  than  1.7  times  the  original  value  for  a2/a  = 


0. 


r»  17  y  »  ■  i  r»  i1  .  ■! 
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Table  1.  Salvo  Kill  Probability 
n  =  2 


\.'v 


ms 


*. .  . 


-v* 

$? 

& 

v 

‘V 


2  .  2 
o^/o. 


o22/a2 

0.00 

0.50 

1.00 

2.00 

4.00 

8.00 

0.0 

1.0000 

0.8333 

0.6667 

0.4667 

0.2889 

0.1634 

0.1 

0.9917 

0.8171 

0.6591 

0.4669 

0.2923 

0.1666 

0.2 

0.9722 

0.7977 

0.6487 

0.4647 

0.2940 

0.1689 

0.3 

0.9467 

0.7767 

0.6365 

0.4609 

0.2946 

0.1706 

0.4 

0.9184 

0.7550 

0.6232 

0.4560 

0.2944 

0.1717 

0.5 

0.8889 

0.7333 

0.6095 

0.4502 

0.2935 

0.1724 

0.6 

0.8594 

0.7120 

0.5956 

0.4439 

0.2920 

0.1728 

0.7 

0.8304 

0.6912 

0.5818 

0.4373 

0.2902 

0.1730 

0 . 8 

0.8025 

0.6712 

0.5681 

0.4305 

0.2881 

0.1729 

0.9 

0.7756 

0.6518 

0.5547 

0.4236 

0.2858 

0.1726 

1.0 

0.7500 

0.6333 

0.5417 

0.4167 

0.2833 

0.1722 
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Table  2  Salvo  Kill  Probability 


n  *  16 


0.00 

0.25 

0.50 

0.75 

1.00 

1.25 

1.50 

1.75 

2.00 

2.25 

2.50 

2.75 

3.00 

3.25 

3.50 

3.75 

4.00 

4.25 

4.50 

4.75 

5.00 

5.25 

5.50 

5.75 
6.00 


0.00  0.50 

1.0000  0.9935 
1.0000  0.9956 
1.0000  0.9965 
1.0000  0.9968 

1.0000  0.9967 
0.9999  0.9962 
0.9997  0.9954 
0.9993  0.9941 

0.9985  0.9922 
0.9972  0.9899 
0.9954  0.9869 
0.9930  0.9833 

0.9900  0.9792 
0.9863  0.9745 
0.9821  0.9692 
0.9772  0.9635 

0.9719  0.9573 
0.9660  0.9507 
0.9597  0.9438 
0.9530  0.9366 

0.9459  0.9291 
0.9386  0.9214 
0.9309  0.9135 
0.9231  0.9055 
0.9151  0.8973 


1.00  2.00 

0.9412  0.7835 
0.9570  0.8234 
0.9659  0.8504 
0.9710  0.8692 

0.9739  0.8824 
0.9752  0.8915 
0.9753  0.8978 
0.9745  0.9017 

0.9728  0.9038 
0.9704  0.9045 
0.9672  0.9038 

0.9635  0.9021 

\ 

0.9591  0.8996 
0.9541  0.8962 
0.9486  0.8921 
0.9427  0.8875 

0.9364  0.8824 
0.9298  0.8769 
0.9228  0.8710 
0.9155  0.8648 

0.9081  0.8583 
0.9004  0.8516 
0.8926  0.8448 
0.8847  0.8378 
0.8767  0.8307 


4.00  8.00 

0.5511  0.3370 
0.6013  0.3784 
0.6399  0.4128 
0.6700  0.4417 

0.6938  0.4662 
0.7126  0.4871 
0.7276  0.5049 
0.7394  0.5201 

0.7485  0.5330 
0.7555  0.5440 
0.7606  0.5534 
0.7642  0.5612 

0.7665  0.5678 
0.7676  0.5732 
0.7677  0.5776 
0.7670  0.5810 

0.7655  0.5837 
0.7634  0.5857 
0.7607  0.5870 
0.7576  0.5878 

0.7541  0.5881 
0.7502  0.5880 
0.7460  0.5874 
0.7416  0.5865 
0.7369  0.5853 


21 
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Experience  shows  that  there  is  one  and  only  one  local  maximum 

2  2  2  2 

of  P  as  a  function  of  o^/a  for  fixed  n  and  a^/a  ,  but  we 
have  not  succeeded  in  obtaining  a  proof. 

The  salvo  kill  probability  can  be  increased  by  increasing 

2  2  2  2 
c^/a  from  zero  for  given  n  and  a^/a  if  the  condition  below 

is  satisfied. 


Proposition  2.1 

The  salvo  kill  probability  P  for  fixed  n,  and  a  is 
not  maximal  at  c?2  =  0  if  and  only  if 


_±_  +  JL_  + 

1+X  2+X  • 


>  1  , 


where 


(2.8) 


.  2.  2 

X  =  a  / a ^ 


Proof 


We  will  show  that 


dP/dp |  >  0 

P=X 


(2.9) 


is  equivalent  to  (2.8).  Using  (2.7),  (2.9)  is  written  as 


n  /  (l-e'"t/A)n”1e"t/A”tdt  <  n  /  (l-e“t//A)  n"1e"t/A~tt/X  dt 


Now  define 


I(x)  =  n  /  (1  -  e”t/A)n"1e~t/X“xt  dt  . 
0 


Then  the  inequality  (2.10)  is  rewritten  as 

1(1)  <  -I'(l)/X  ,  (2.1 


where  I  * (1)  means  the  value  of  the  derivative  dl/dx  at  x  = 
Let  s  =  exp(-t/X)  in  the  integral,  then  I (x)  is 


I(x)  =  nX  f1  sAx(l-s)n“1  ds 
0 

,  n  n-1  1 

”  Xx+1  Xx+2  ’ ‘ ’  Xx+n  * 


Therefore, 


I' (x)  _  _  r  X 

l(x)  “  ^1  Xx+3 

The  inequality  (2.8)  is  readily  obtained  from  this 
equation  with  x  =  1  and  (2.11) .  □ 

Proposition  2 . 2 

A  sufficient  condition  that  the  salvo  kill  probability 
P  for  fixed  n,  o. ,  and  a  is  not  maximal  at  a,  =  0  is 


n/(e-l)  -  1  . 


(2.12) 


2,2  < 

a  /a1  = 


Proof 

Since  l/(x+A)  is  decreasing  in  x. 


n  1 

,L  & 


> 


n 


3+1 

j 


dx 

x+A 


n+1  , 

f  dx 

J  ^  x+A 


In 


n+l+A 

1+A 


If  the  right-hand-side  is  equal  to  or  greater  than  unity, 
i.e.,  l+n/(l+A)  =  e,  then  the  condition  (2.8)  in  Proposition 
2.1  is  satisfied.  The  relation  (2.12)  is  immediate  from  the 
above-mentioned  inequality.  □ 

Proposition  2.3 

Assume  that  (2.8)  holds.  A  necessary  condition  that  the 
salvo  kill  probability  be  maximum  at  some  point  a2  >  0  for 
fixed  n,  a1  and  a  is 


?  (-1)  j“1(n)j(^) j"1  P+3 7.1.  =  0  ,  (2.13) 

j=l  3  p  (p+j) 


where 


A  =  a 2/al  ,  p  =  (c i2+o22)/o\  . 


Proof 

Differentiate  (2.7)  with  respect  to  p,  and  equat  it  to 
0.  The  equation  (2.13)  is  immediate.  □ 


Example.  Case  n  =  2 


Equation  (2.13)  is  generally  too  complicated  to  be  solved 
for  a 2/ ot  in  closed  form.  Here  we  solve  it  for  the  simplest 
case,  n  =  2. 

The  inequality  (2.10)  with  n  =  2  is 


1 

1+A 


+ 


1 

2+A 


> 


1  , 


and  after  a  bit  of  algebra,  this  reduces  to 

1/A  =  aj/a2  >  2/(/5  -  1)  .  (2.14) 

2  2 

Now,  assume  that  a-^/a  is  large  and  satisfies  (2.14).  From 
(2.13) 


2  _ £ _  -  2  —  _  =  0 

(p+D2  p  (P+2r 


Therefore, 


2/ 

<W  a 


2 


p/A  -  1 


(p+1)  3 
p (p+2) 2 


1  , 


Oi/a 


2 


1/A 


(P+1)3 

- T  • 

P  (p+2)^ 


(2.15) 


2  2 

This  is  a  parametric  expression  of  maximizing  c^/oi  as  a 
2  2 

function  of  o^/ol  the  curve  of  which  is  shown  in  Fig.  2.  The 
value  of  the  parameter  p  which  corresponds  to  the  point 


(a^/a  =  (/5’  +  l)/2,  a^/ a  =0)  is  p  =  pQ,  where 

pQ  =  </5  -  l)/2  . 

When  p  approaches  zero, 

oj/a2  ~  1/4 p 2  ,  a2/a2  -  1/4 p  , 

and  therefore,  the  asymptotic  form  of  the  curve  is 

o2/a2  *  j  a-j/ct  .  (2.16) 

B.  COMPUTING  THE  SALVO  KILL  PROBABILITY 

In  the  preceding  section,  the  Weiss  formula  for  the  salvo 
kill  probability  was  given  in  (2.7) .  According  to  Breaux 
[Ref.  7],  this  formula  is  not  suitable  for  calculation  for  a 
large  n  because  it  is  an  alternating  series.  Breaux  pointed 
out  that  the  salvo  kill  probability  could  be  expressed  in 
terms  of  an  Incomplete  Beta  Function,  and  that  for  large  n, 
Tang's  method  [Ref.  8]  would  be  advisable  for  calculating  an 
Incomplete  Beta  Function. 

First,  we  present  a  recurrence  formula  for  calculating  the 
salvo  kill  probability  which  is,  in  principle,  the  same  as 
that  of  Tang's  method.  Its  generalization  will  be  given  in 
a  later  chapter. 


Proposition  2.4 

Let  us  denote  the  salvo  miss  probability  by  n  weapons  as 
Q(n)  =  1-P . 


Q(n)  =  /  (1  -  ye"t/p)n  e_t  dt  . 

0 


(2.17) 


Then  Q(n)  satisfies  the  following  recurrence  equations: 


Q  (n)  = 


_S_Q(n-i)  +  , 


Q (0)  =  1  , 


n  =  1,2, 


(2.18) 


(2.19) 


where 


a2/  (a2  +  a^)  / 


Ca2  +a2)/a2 


Proof 


From  (2.17) 


Q(n)  =  /  (1  -  ye“t/p)n  e~t  dt 

0 


=  /  (l-ye-^jH-V^t-u/  a-we‘t/p)n"1«“‘t/p‘tdt 
0  0 
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The  first  term  in  the  right  hand  side  is  Q(n-l).  In  the 
second  term,  we  integrate  it  by  parts,  and  have 

Q  (n)  *  Q(n-l)  -  £(l-ye“t/p)  ne_t  |  -  £  /  (l-ye_t/p)  ne-t  dt 

n  0  n  0 

=  Q(n-l)  +  £<l-y)n  -  £•  Q  (n)  . 

The  equation  (2.18)  is  immediate.  Q 

It  is  worthwhile  to  note  that  the  method  presented  in 
Proposition  2.4  gives  not  only  the  value  of  Q(n)  for  a 
specified  nQ,  but  also  gives  all  the  Q(n) 's  up  to  nQ.  It 
is  an  important  advantage  of  this  method  over  the  Weiss 
formula  because  we  are  often  interested  in  the  salvo  kill 
probabilities  for  several  n's.  As  is  shown  in  Table  3,  the 
computational  time  for  Q(n)  by  our  recurrence  algorithm  is 
about  twice  the  time  required  for  computation  of  a  single 
Q(n)  using  the  Weiss  formula.  When  the  salvo  kill  probabili¬ 
ties  for  more  than  a  single  n  are  needed,  therefore,  use  of 
the  recurrence  algorithm  given  in  Proposition  2.4  is 
recommended . 

We  have  another  efficient  algorithm  for  computing  the 
salvo  kill  probability,  which  is  good  even  for  large  n. 
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Table  3.  Comparison  of  the  Processing  Time  for 
Computing  the  Salvo  Kill  Probability 


Mean 

time  (microsecond) , 

p  =  1,  and  y 

• 

m 

• 

o 

II 

n 

Weiss  formula 

Recurrence 

formula 

10 

92 

175 

20 

189 

383 

30 

300 

587 

Mean 

time  (microsecond) , 

p  =  1,  and  n 

=  10. 

V 

Weiss  formula 

Recurrence 

formula 

0.1 

90 

196 

0.5 

92 

175 

0.9 

100 

186 

Computation  is  on  the  IBM  3033  installed  in  the  Naval 
Postgraduate  School.  Each  value  is  the  mean  of  4  trials. 
One  trial  is  1000  repetitions  of  calculation,  and  the 
elapsed  time  units  are  divided  by  1000  and  multiplied  by 
26  (microseconds) . 


Proposition  2.5 


i 


The  salvo  miss  probability  Q(n)  by  n  weapons 


Q (n)  =  /  (1  -  ye"t/p)n  e-t  dt 

0 


(2.17) 


satisfies  the  following  relationship 


Q(n)  =  l  (?)  (l-y)n_k  yk  q(k)  , 


(2.20) 


where  q(k)  is  the  salvo  miss  probability  by  k  weapons  with 
y  =  1/ 


oo 

q  (k)  =  /  (1 -e't/p)k  e”*  dt 

0 


Proof 


Q(n) 


/  { (1  -  y)  +  y(l-e"t/p)}n  e”fc  dt 
0 


¥  ,n>  .  n-k  k  f  ..  -t/p^K  -t  .. 

I  (v)  (1-y)  y  /  (1  -  ye  '  )  e  dt 

k=0  *  0 


Corollary  2.6 


The  salvo  miss  probability  (2.17)  can  be  calculated  by 
the  formula 


n 

Q(n)  =  l  A.  ,  (2.21) 

k=0 


where 


Ak^Ak-l  =  y (n-k+1) /(1-y) (k+p)  ,  k  =  1,2,... 
A0  =  (l-y)n  . 


(2.22) 


Proof 

As  can  be  seen  from  (2.18)  and  (2.19)  with  y  =  1,  q(0) 
is  unity,  and  q(k)/q(k-l)  is  k/(k+p).  Utilizing  the  fact 
that  (J)/^)  =  (n-k+l)/k,  the  corollary  can  be  obtained  by 
substitution.  □ 

Note  that  all  the  terms  in  the  expression  (2.24)  are 
positive,  so  there  is  no  problem  arising  from  cancellation 
of  terms  with  alternating  signs  as  was  observed  in  the  Weiss 
formula  (2.7)  . 

In  checking  the  computation,  the  following  P's  for  special 
cases  would  be  of  use. 

(1)  The  case  =  0  corresponds  to  p  -*•*>,  and  we  have 


(2)  The  case  =  0  corresponds  to  y  =  1,  p  =  A,  and  so 
from  the  recurrence  formula,  we  have 


n(n-l)  1 _ 

(n+A) (n-l+A)  •••  (1+A) 


(2.24) 


2  2  2 

(3)  The  case  a  +  a2  =  ai*  The  formula  (2.5)  with  p  =  1 
is  easily  integrated  to  give 

P  =  1  -  -  (1-A)n+1)  ,  (2.25) 

where 

U  =  a2/ (a2  +  o\)  ,  A  =  a 2 /a2  . 

C.  BOUNDS  TO  THE  SALVO  KILL  PROBABILITY 

In  Chapter  I  it  was  pointed  out  that  an  upper  bound  of 
the  salvo  kill  probability  is  given  by  the  kill  probability 
of  n  repeated  independent  shots. 

P  <  P 


where 

00  00 

P  ■  1  -  [/  /  (1  -  p(u,v) }f.  (u,v) dudv]n  . 

—00  —  00 

Under  the  assumptions  of  the  present  chapter,  we  have 
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Proposition  2.7 

An  upper  bound  to  the  salvo  kill  probability  is  given  by 

P  =  1  -  (1  -  ^ )n  ,  (2.26) 

where 

2  2  2  2  2 
X  =  a  ,  p  =  (a  +0^)70^  . 

As  to  a  lower  bound,  we  have 

Proposition  2.8  (Merritt  and  King) 

A  lower  bound  to  the  salvo  kill  probability  is  given  by 

P  *  1  -  p(nu)"p  rny(p)  ,  (2.27) 

where  r  (p)  is  the  Incomplete  Gamma  Function, 

X  - 

r  (p)  =  /  e~C  tp“A  dt  , 

x  0 

and 

y  =  a2/ (a 2  +  o\)  ,  p  =  (a 2+a2)/a2  . 

Table  4  and  Table  5  show  how  the  true  kill  probabilities 

are  bracketed  by  these  two  bounds.  The  triplets  of  entries 
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Table  4  Salvo  Kill  Probability  and  Its  Bounds 

n  =  2 


& 


2  .  2 
aL/a 


2.  2 
a0/a 

0 . 00 

0.50 

1.00 

2.00 

4.00 

3.00 

0.8889 

0.7500 

0.6400 

0.4898 

0.3306 

0.1994 

0.5 

0.8889 

0.7333 

0.6095 

0.4502 

0.2935 

0.1724 

0.7364 

0.6187 

0.5216 

0.3911 

0.2580 

0.1527 

0.7500 

0.6400 

0.5556 

0.4375 

0.3056 

0.1900 

1.0 

0.7500 

0.6333 

0.5417 

0.4167 

0.2833 

0.1722 

0.6321 

0.5443 

0.4715 

0.3679 

0.2532 

0.1552 

0.6400 

0.5556 

0.4898 

0.3951 

0.2840 

0.1814 

1.5 

0.6400 

0.5524 

0.4325 

0.3829 

0.2696 

0.1689 

0.5507 

0.4832 

0.4267 

0.3429 

0.2442 

0.1542 

0.5556 

0.4898 

0.4375 

0.3600 

0.2653 

0.1736 

2.0 

0.5556 

0.4881 

0.4333 

0.3524 

0.2554 

0.1643 

0.4866 

0.4333 

0.3882 

0.3192 

0.2338 

0.1516 

0.4898 

0.4375 

0.3951 

0.3306 

0.2489 

0.1664 

2.5 

0.4898 

0.4365 

0.3925 

0.3255 

0.2418 

0.1593 

0.4353 

0.3922 

0.3554 

0.2977 

0.2233 

0.1482 

0.4375 

0.3951 

0.3600 

0.3056 

0.2344 

0.1597 

3.0 

0.4375 

0.3944 

0.3583 

0.3021 

0.2292 

0.1542 

0.3935 

0.3580 

0.3274 

0.2784 

0.2131 

0.1444 

0.3951 

0.3600 

0.3306 

0.2840 

0.2215 

0.1536 

3.5 

0.3951 

0.3596 

0.3294 

0.2815 

0.2175 

0.1492 

0.3588 

0.3291 

0.3032 

0.2611 

0.2034 

0.1404 

0.3600 

0.3306 

0.3056 

0.2653 

0.2099 

0.1479 

4.0 

0.3600 

0.3303 

0.3048 

0.2635 

0.2068 

0.1443 

0.3297 

0.3045 

0.2823 

0.2457 

0.1944 

0.1365 

X3, 

jsS 


Triplets  are,  from  the  top,  P,  the  upper  bound, 
P,  the  true  value,  and  P,  the  lower  bound. 
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Table  5  Salvo  Kill  Probability  and  Its  Bounds 

n  *  16 


a* /cl 


2 


0.00 

0.50 

1.00 

1 . 0000 

0.9997 

0.9985 

1.0000 

0.9967 

0.9739 

0.9997 

0.9944 

0.9688 

0.9985 

0.9954 

0.9900 

0.9885 

0.9922 

0.9728 

0.9952 

0.9862 

0.9644 

0.9900 

0.9821 

0.9719 

0.9900 

0.9792 

0.9591 

0.9817 

0.9685 

0 . 9469 

0.9719 

0.9597 

0.9459 

0.9719 

0.9573 

0.9364 

0.9592 

0.9432 

0.9215 

0.9459 

0.9309 

0.9151 

0.9459 

0.9291 

0.9081 

0.9305 

0.9130 

0.8917 

0.9151 

0.8987 

0.8819 

0.9151 

0.8973 

0.8767 

0.8983 

0.8802 

0.8597 

0.8819 

0.8650 

0.8481 

0.8819 

0.8639 

0.8441 

0.8647 

0.8468 

0.8273 

0.8481 

0.8313 

0.8147 

0.8481 

0.8305 

0.8117 

0.8310 

0.8136 

0.7953 

2.00 

4.00 

8.00 

0.9900 

0.9459 

0.8147 

0.8824 

0.6938 

0.4662 

0.8750 

0.6867 

0.4611 

0.9719 

0.9151 

0.7824 

0.9038 

0.7485 

0.5330 

0.8935 

0.7388 

0.5258 

0.9459 

0.8819 

0.7515 

0.8996 

0.7665 

0.5678 

0.8864 

0.7546 

0.5590 

0.9151 

0.8481 

0.7222 

0.8824 

0.7655 

0.5837 

0.8674 

0.7522 

0.5739 

0.8819 

0.8147 

0.6945 

0.8583 

0.7541 

0.5881 

0.8424 

0.7402 

0.5778 

0.8481 

0.7824 

0.6684 

0.8307 

0.7369 

0.5853 

0.8145 

0.7229 

0.5749 

0.8147 

0.7515 

0.6439 

0.8016 

0.7167 

0.5780 

0.7857 

0.7030 

0.5677 

0.7824 

0.7222 

0.6209 

0.7724 

0.6952 

0.5680 

0.7570 

0.6819 

0.5579 

Triplets  are,  from  the  top,  P,  the  upper  bound, 
P,  the  true  value,  and  P,  the  lower  bound. 
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in  the  table  are,  from  the  above,  the  upper  bound  P,  true 

value  P,  and  P,  the  lower  bound.  It  is  observed  that 
_  2  2 

P-P  is  not  small  for  smaller  values  of  c^/a  ,  but  gets  smaller 

2  2.  2  2  — 
when  / a  is  increased.  If  a^/a  is  small,  p  and  P  are 

close.  On  the  other  hand,  P  is  fairly  close  to  P  over  wide 
2  2  2  2 

ranges  of  c^/a  and  c^/cx  .  Figures  3  and  4  illustrate  this. 

It  seems  reasonable  that  Merritt  and  King  used  P  as  an 
approximate  formula  for  p  in  their  study  on  the  optimal 
ballistic  dispersion,  because  computation  was  really  a 
problem  at  that  time.  Nowadays,  however,  direct  computation 
of  P  is  much  more  natural. 


III.  OPTIMIZATION  OF  SALVO  FIRING 


In  Chapter  II  we  have  learned  the  following:  The  salvo 
kill  probability  is  a  decreasing  function  of  a  if  other 
parameters  are  kept  constant.  Therefore,  to  obtain  higher 
salvo  kill  probability,  it  is  necessary  to  make  as  small 
as  possible.  The  salvo  kill  probability  as  a  function  of 
a2  is,  however,  complicated.  If  c-^  is  sufficiently  small, 
the  salvo  kill  probability  is  monotone  decreasing  in  o2, 
but  if  a1  is  not  so  small  and  satisfies  the  condition  (2.8) , 
then  the  salvo  kill  probability  is  not  maximal  at  a2  =  0, 
but  takes  a  maximum  value  at  some  a2  >  0.  In  other  words,  the 
smallest  ballistic  dispersion  is  not  always  most  desirable. 
Deliberate  increase  of  the  ballistic  dispersion  is  encoun¬ 
tered  in  shotguns  and  in  related  military  weapons.  In  this 
chapter  we  present  an  approximate  formula  for  the  optimal 
ballistic  dispersion.  We  will  first  briefly  sketch  Walsh's 
theory  [Ref.  9]  since  our  approximate  formula  is  based  on 
it.  In  Section  B  our  approximate  formula  is  presented. 

Its  accuracy  will  be  studied  in  Chapter  IV. 

A.  WALSH  THEORY 

The  Walsh  model  is  based  on  assumptions  similar  to  those 
of  Chapter  I.  A  salvo  of  n  weapons  is  fired  at  a  target. 

The  impact  point  error  is  composed  of  two  parts,  a  random 
bias  which  is  common  to  all  n  weapons,  and  a  round-to- round 
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error.  A  Cartesian  coordinate  system  is  chosen  so  that  its 
origin  coincides  with  the  center  of  impact  of  n  weapons. 

The  coordinates  of  the  random  position  of  the  target  are 
denoted  as  U  and  V  with  respect  to  this  system.  The  joint 
probability  density  of  U  and  V  is  given  by  f^(u,v) .  The 
impact  points  of  weapons  (X^,Y^)  are  independent  and  iden¬ 
tically  distributed  random  variables  with  joint  density 
f2 (x,y) .  It  is  assumed  that  impact  of  a  weapon  at  (xry) , 
given  the  target  is  at  (u,v) ,  destroys  the  target  with 
probability  D(u-x,v-y)  independently  of  other  weapons. 

The  lethal  area  of  a  weapon  is  denoted  as  A: 


oo  OO 

J  j  D(x,y)dxdy  «  A  . 

— OO  —00 


The  salvo  kill  probability  is  given  by  (1.2), 


P  =  1  -  /  /  {1-  p(u,v) jHf^ (u,v)dudv  ,  (1.2) 

—00  —00 


where  p(u,v)  is  the  conditional  destruction  probability  by 
a  single  weapon  given  that  the  target  is  at  (u,v) . 


00  00 

p (u, v)  *  /  /  D(u-x,v-y) f2 (x,y) dxdy  .  (1.1) 

—00  —00 


Note  that  p(u,v)  can  be  integrated  to  give  A. 


A  . 


(3.1) 


/  /  p (u,v) dudv  =  A  /  /  f  (x,y)dxdy  = 


—CO  —00 


—00  —00 


The  problem  is  to  derive  the  optimal  ballistic  disper- 
2 

sion,  a 2-  Instead  of  solving  the  original  problem,  Walsh 

solved  a  revised  problem  that  was  different  from  the  original 

2 

one  in  that  the  variable  is  not  c2  itself  but  rather  the 
function  p(u,v) . 

Proposition  3.1  (Walsh) 


max  P 
p(  •  • ) 


uu  uu 

■  1  -  /  /  { 1  -  p(u,v) }nf ^ (u, v) dudv  (3.2) 


subject  to 


/  /  p(u,v)dudv  ■  A  , 


—00  —00 


p  (u,  v)  =  0 


has  a  solution 


(3.3) 


p* (u, v)  = 


n  _  r  C  i  n-J 
1 


(u,v)  e  E 


Cu,v)  e  E 


C3.41 


where 


42 


E 


{ (u,v) 


(u,V) 


> 


C) 


( 


C  is  positive  and  is  given  by 


/  /  p*(u,v)  dudv  =  A  .  ( 

E 


The  corresponding  salvo  kill  probability  is 


P*  =  CA  +  /  /  {f,  (U/V)  -  C}  dudv  .  ( 

E 

In  case  the  random  target  bias  is  given  by  a  normal 
distribution,  we  also  have 

Corollary  3.2 

00  00 

max  P  *  1  -  /  /  {1-  p  (u,v)  }nf  ^  (u,v)  dudv  (■ 

p  (  •  •  )  — 00  — CO 

subject  to 


ao  oo 

/  /  p(u,v)  dudv  =  A  , 

—00  —00 

( 

p(u,v)  >  0  , 


where 


2ir  a  o 


( 


has  the  solution 


F-xv: 


v.v. 


h>; 


K-; 


B 


N>"v' 

W 


v'V 

k^v  ' 


Iftgj 

p.»?>r 


1  2  2 

1  -  exp{-<j>  +  2Xn^Ty(~T  +  (u,  v)  e  E 


p*(u,v) 


a  a 

U  V 


(3.11) 


(u,v)  e  E 


where 


E  «  {(u,v)|u2/o2  +  v2/a2  i  2(n-l)<t>) 


4>  is  positive  and  is  determined  by 


-4> 


e  T  -  1  +  <j)  -  a/2tt  (n-1)  °u  =0 


(3.12) 


The  maximum  value  of  the  salvo  kill  probability  is  given  by 


p*  =  i  _  e-(n-l)  <f>{1  +  (n-1)  Cl  -e"^)} 


(3.13) 


Proof 

The  equation  (3.11)  is  readily  obtained  from  (3.4)  in 
Proposition  3.1  with  f1(u,v)  given  by  (3.10)  and  C, 


C  = 


-(n-1)  <j> 


2*  »u  ®v 


(3.14) 


Substituting  (3.11)  in  (3.6)  and  changing  variables  into 
(r, 0}  by 
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'  .  V  ’ 

t  VI 


u  =  /2  (n-1)  au  r  cos  9  ,  v  =  /2  (n-1)  a^r  sin  9  (3.15) 

we  have 


A  =  /  /  p*(u,v)  dudv 

E 

A  ,  _2 

=  2tt  (n-1)  ouav  /  (1  -  e”^  )  2r  dr 

*  2ir  (n-1)  au  av  ( <J>  -  1  +  e-^)  . 


Substituting  (3.10)  and  (3.14)  in  (3.7)  with  (3.15) , 


j*  = 


A  e~  (n-1)  <p  +  2tt 


2ir a  a 
u  v 


so:*”-11  I*  (e-'"-1'  r2-e-  (n-1)  * 


u  V 


=  1  -  +  (n-1)  4>  -  A/2TTU  a  }  . 


U  V 


This  formula  together  with  (3.12)  gives  (3.13).  □ 

For  simplicity,  let  us  call  the  function  p*(u,v).  the 
(Walsh)  optimal  p  function,  and  the  associated  salvo  kill 
probability  P*  the  Walsh  optimal  P.  It  should  be  noted  that 
the  constraint  in  the  revised  problem  is  weaker  than  that  of 
the  original  problem,  and  so  the  Walsh  optimal  P  is  an  upper 
bound  to  the  maximum  salvo  kill  probability  of  the  original 
problem.  It  should  also  be  noted  that  the  Walsh  optimal  p 
function  is  not  feasible  in  the  original  problem  (being  non 
zero  on  only  a  finite  set)  if  dispersion  is  present. 


)  2rdr 
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B.  APPROXIMATE  FORMULA  FOR  THE  OPTIMAL  BALLISTIC  DISPERSION 
Walsh  developed  the  theory  stated  above,  but  he  did  not 
extend  his  theory  beyond  the  results  (3.4) -(3. 7).  The  impor^ 
tant  part  left  by  Walsh  and  subsequent  workers  is  to  bridge 
the  gap  between  the  revised  problem  and  the  original  prob¬ 
lem.  Walsh  got  the  optimal  p  function  as  given  in  (3.4), 
whereas  we  need  the  optimal  density  function  f2(x,y).  The 
material  in  this  section  bridges  that  gap. 

Let  us  define  three  double  Fourier  transforms  as  follows; 

~  00  oo  i  ((Jj^u-hu  „v) 

P*  ((iK/OJj)  ■  /  /  e  p*(u,v)dudv, 

— oo  — oo 

~  oo  oo  i  (oj  x-Hij^y) 

D(u1,u2)  *  /  /  e  1  Z  D(x,y)dxdy,  (3.16) 

— OO  —00 

oo  oo  i(u>.x+w2y) 

f2(oj1,w2)  -  /  /  e  f 2  (x,y)  dx  dy  . 

Then  from  (1.1)  we  have 

p*(u1,0J2)  *  D(w1,ta)2)  f2  (w^,w2)  .  (3.17) 

In  the  Walsh  theory,  D(x,y)  is  a  given  function,  and 
the  Walsh  optimal  p  function  is  given  by  Proposition  3.1. 
Therefore,  both  D(oj^,w2)  and  p*  (uj^ ,u»2)  in  (3.17)  are  now 
known  functions.  Thus  the  Fourier  transform  of  the  "optimal" 
f~(x,y)  can  be  calculated  by 


P*(CJ,  ,w9) 

f2(w1#w2)  =  r - - — —  ,  (3.18) 

D((jl)1,U2) 

* 

and  we  will  get  f2(x,y)  as  its  inverse  Fourier  transform. 

This  might  have  been  the  end  of  Walsh  theory,  but  there  are 
diff iculties.  First,  we  do  not  have  any  guarantee  that 
f2(x,y)  is  a  density  function,  because  Walsh  solved  the 
restricted  optimization  problem  subject  to  p(u,v)  >  0  but  not 
f2(x,y)  >  0.  Further,  even  if  this  problem  were  solved,  there 
would  still  remain  the  problem  of  feasibility  of  the  f2  func¬ 
tion.  Walsh  suggests  in  his  paper  that  it  would  be  possible 
to  obtain  a  real-world  dispersion  function  f2(x,y)  as  the 
theory  requires,  at  least  approximately,  by  careful  design 
of  the  ammunition,  but  this  may  be  asking  too  much  of  the 

manufacturing  process.  We  assume  below  that  the  function 

2 

f2(x,y)  is  normal,  with  only  the  variance  a 2  subject  to 
manufacturing  control . 

The  problem,  then,  is  to  determine  or  approximate  the 

2 

optimal  ballistic  dispersion  a2,  given  the  damage  function 
D (x,y)  and  the  function  p*(x,y)  from  (say)  Proposition  3.1. 

In  the  following,  both  D(x,y)  and  p*(x,y)  are  assumed  to  be 
even  functions  of  x  and  y  respectively,  and  to  have  "moments" 
of  all  degrees.  The  lethal  area  is  denoted  as  A. 

00  00  OO  00 

/  /  xD(x,y)dxdy  =  /  /  yD(x,y)dxdy  =  •••  =  0  ,  (3.19) 

— OO  —00  —00  —00 
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«  0  , 


/  /  xp*(x,y)dxdy  =  j  j  yp*(x,y)dxdy  =  ••• 

—00  —CO  —00  — oo 


00  00 

/  /  D(x,y)  dxdy  =  A  , 

— oo  — oo 


/  /  x2kD(x,y)  dxdy  <  ®,  k  =  1,2,  ... 

—CO  —co 


/  /  yl  p*(x,y)  dxdy  <  <*> ,  k 


1,2, 


(3.20) 


etc. 

For  the  sake  of  generality,  we  actually  consider  the  non 
symmetric  case  where  the  optimal  dispersion  may  be  differ¬ 
ent  in  the  x  and  y  directions.  Let  the  two  dispersions  be 
2  2 

a  and  a  .  Then 
x  y 

f2<x’y>  ■  jwV  e 

x  y 

with  the  double  Fourier  transform 


f2(a)1,a)2)  = 


=  e 


2  2  2  2 
xul  °y“2 

~~i - V 


C  3 . 22). 


Expanding  the  exponential  function  in  the  double  Fourier 


transform  of  D(x,y) ,  and  integrating  term  by  term,  we  have 


00  00  i(w,x+u>2y) 

D(oj^,cj2)  =  /  /  e  D(x,y)dxdy 


— oo  — oo 


OO  ,  j  CO  CO 

l  XT  I  /  (uj,  x+ui-y)  (x,y)  dxdy 
j»0  3*  -  “  1  1 


—00  —00 


10,  ,  »  00  9 
1  1  (  <  „2, 


=  A[1  -  -2~  i  /  /  x  D(x,y)dxdy 

—03  —00 

2 

(o9  ,  ®  “>  9 

- 5~  r  /  /  y  u  (x,y)  dxdy  +  ...] 


—00  —  oo 


where  the  third  line  is  obtained  because  D  is  even.  There- 

2  2 

fore,  for  sufficiently  small  values  of  to1  and  u>2. 


D(uj^,u)2)  =  Dq(o}^,u)2) 


(3.23) 


where 


Dq  (oj1,u)2) 


A  exp  [-- y  i  /  /  x2D (x,y)  dxdy 


*.00  —00 


2 

-  x  /  /  y2D(x,y)  dxdy] 


—00  —00 


(3.24) 


Similarly, 
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(3.25) 


p*(w1,a)2)  =  pQ(ui1,u2) 


2  2 

for  small  and  uj2,  where 


W1  1  00  °° 


Po(a)l'a>2)  *  A  exP  I - 2"  K  f  /  x  p*(x,y)  dxdy 


—  CO  — oo 


-  T  f  /  Y^P*(x,y)  dxdy]  . 


2  A 


(3.26) 


—00  —00 


2  2 

Now  we  determine  two  parameters  and  in  (3.21)  by 
the  relation 


Pn(w1»u?) 

f2  (u.  ,cO  =  J 


(3.27) 


Then  from  (3.22),  (3.24),  and  (3.26)  we  get  an  approximate 
formula  for  the  optimal  ballistic  dispersion: 


O  oo  oo 

°x  *  /  /  x2{p* (x,y) -D(x,y)  }  dxdy,  0]  + 


—00  —00 


(3.28) 


00  00 


oj  -  [jj-  /  /  y2{p*(x,y) -D(x,y) }  dxdy,  0]  + 

— O 


-00  —00 


where 


[a, 0  ] + 


a  if  a  *  0 
0  otherwise 


(3.29) 
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Note  1 


In  this  proposition,  the  approximate  formula  for  the 

optimal  dispersion  is  derived  from  the  condition  that 

f2(ajl/ij2^  coincides  with  f*(u  tu>2)  for  infinitesimally  small 
2  2 

and  u^.  It  would  be  better  if  we  could  derive  the  approxi 
mate  formula  from  the  condition  that  f 2^1 '“2)  is  neariY 
equal  to  f2  (00^, u^)  over  a  wider  domain  of  and  which 
include  (0,0)  .  In  the  next  chapter  we  will  give  an  idea  of 
improving  the  approximate  formula  taking  this  into  account. 

Note  2 

Generally,  approximation  (3.23)  with  (3.24)  is  accurate 
only  up  to  the  second  power  of  and  u^.  If  the  damage 
function  is  of  Gaussian  type, 


then 


D0(w1,w2)  =  D(u1,u2) 

is  an  exact  relationship. 

As  a  corollary,  we  have  the  following  approximate  formula 
for  the  model  with  normal  errors  and  Gaussian  damage  function, 

Consider  a  salvo  firing  of  n  weapons.  The  density  func¬ 
tions  of  the  random  bias  and  the  round-to-round  ballistic 


error  are,  respectively. 


£1(u,v)  = 


2tt  a  a 
u  v 


e 


2a  2a' 
u  T 


and 


f2(x,y)  = 


2tt  a  a 
x  y 


2 

2a“  2a2 

x  y 


the  damage  function  is  given  by 


D (u-x, v-y)  =  e 


(u-x) 2  (v-y) 2 

- 5 - ? — 

2*u  2ov 


(3.32) 


(3.33) 


(3.34) 


Then  approximate  formulae  for  the  optimal  ballistic  dispersion 
for  fixed  n,  au,  av,  ax,  and  ay  are 

(ax/0tx)2  =  [l(n-1)2^(Vax)3(Vay)“(n‘1)(ou/ax)2"1'°]  + 

(3.351 

<VV2  =  fl(n-l)24>2(VV3(au/ctx)-(n'1)  (Vay)2_l'0J  + 

where  <p  is  determined  by  the  equation 

e"*  -  1  +  <0  -  (ax/au)  (ay/av)/(n-l)  =  0  .  (3.36) 
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An  upper  bound  to  the  maximum  salvo  kill  probability  is  given 
by 


P  =  1  -  e"(n“1)4){l  +  (n-1)  (1  -  e~*) }  .  (3.37) 

The  solution  to  the  Walsh  problem  with  the  above  men¬ 
tioned  assumptions  is  already  given  in  Corollary  3.2.  By 
use  of  (3.11),  we  have 


00  oo 


jr  /  /  u2p*(u,v)  dudv 


—CO  — oo 


o  o  /2(n-l)$  ,  0  -<p  +  4  v 

u  v  2it  /  4 r2  {1  -  e  2(n'1)}  rdr 


u  V, 


2.2 


A  -{ (n-1)  <p  -  (n-l)A/Ttauav} 


and 


00  oo 


r  f  I  x2D(x,y)  dxdy  = 


—00  —00 


2 

ax  * 


Substituting  these  relations  into  (3.28),  we  get  the  approxi 

2 

mate  formula  for  ox  (3.35) .  Similar  arguments  follow  for 


SIMPLE  SALVO  MODEL  REVISITED 


In  this  chapter  we  apply  the  general  approximate  formula 

for  the  optimal  ballistic  dispersion  obtained  in  the  last 

chapter  to  the  Simple  Salvo  Model  dealt  with  in  Chapter  II 

and  investigate  its  accuracy. 

The  assumptions  in  this  chapter  are  the  same  as  those  in 

Chapter  II.  We  consider  a  salvo  of  n  weapons  against  a 

target.  The  impact  point  error  is  composed  of  two  parts, 

the  random  bias  and  the  round-to-round  error.  The  position 

of  the  target  (U,V)  has  a  circular  normal  distribution  centered 

2 

at  (0,0)  with  variance  o1 .  The  impact  points  of  the  ith 
weapon  (X,Y)  are  assumed  independent  and  identically  dis¬ 
tributed  circular  normal  random  variables  centered  at  (0,0) 
with  variance  a^.  Finally,  we  assume  the  circular  Gaussian 
damage  function  with  parameter  : 


2 ,  2 
u  +v 


1  1 
f.  (u,v)  =  y  e 

1  2-na* 

.  si±zi 

i  2a\ 

f2(x,y)  =  5-  e 

2* 


(4.1) 


(4.2) 


(u-x) 2+(v-y) 2 


D (u-x, v-y)  =  e 


(4.3) 


In  (3. 35) ,  let 


a 


u 


a. 


a 


(4.4) 


then  we  have  the  following: 

An  approximate  formula  for  the  optimal  ballistic  disper¬ 
sion  of  the  Simple  Salvo  Model  is  given  by 


2.  2  , .2  2  .  nl+ 
02/a  =  I<f>  z  -  z-1,0] 


(4.5) 


where 

z  =  (n-1)  a^/a2  *  l/Ce"*  -  1  ♦  $)  .  (4.6) 

An  upper  bound  to  the  maximum  salvo  kill  probability  is  given 
by 


P  *  1  -  e_(n"1)4){l  +  (n-1)  (1 -e"*)  }  .  (4.7) 

2  2  2  2 
The  true  optimal  o2/a  is  a  function  of  n  and  a^/a  ,  but 

2  2 

according  to  (4.5)  and  (4.6)  the  approximate  a2/a  is  deter- 

2  2 

mined  by  z  =  (n-ljo^/o  ,  $  being  a  parameter.  To  obtain  an 

2  2 

explicit  formula  for  o^/ci  in  z  instead  of  the  parametric 
form,  we  have 
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Lemma  4 . 1 


The  series  solution  to  the  equation  in  <0, 
e_<J>  -  1  +  <0  -  u2/2  =  0 


(4.8) 


is  given  by 


u2  +  u3/3  +  u4/12  +  u5/60  + 


(4.9) 


Proof 

Let  the  left-hand-side  of  (4.8)  be  expanded  in  <j>. 

Then 


u2  =  4>2(1  -  <0/3  +  <o2/12  -  <03/6O  +  ...)  . 


If  the  terms  in  the  right-hand- side  of  (4.9)  are  expanded 
in  <0,  we  have 


2 


u 


u4/12 

u5/60 


<02  -  <03/3  +  <04/12  -  <05/6O  +  ... 
<03/3  -  <04/6  +  <05/18  -  ... 

<04/12  -  <05/18  +  ... 

<05/6O  -  . .  . 


from  which  (4.9)  is  obvious. 


□ 
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Using  the  Lemma  and  the  relationship  (4.6),  we  can  express 
2  2 

02/01  in  (4.5)  as  a  power  series  in  /2/z .  If  we  keep  only 
the  first  three  terms  in  the  expansion,  we  get  an  alterna¬ 
tive  approximate  formula  for  the  optimal  ballistic  disper- 

2 

sion  which  we  will  denote  as  a  . 

C l 

a2/a2  =  |  /572  -  f  +  3^-  /375T  ,  (4-10) 


with 


z  =  (n-1)  a2/ct2  . 


In  Table  6  two  values  of  the  approximate  optimal  ballis¬ 
tic  dispersion  are  compared.  One  is  calculated  with  (4.5) 
and  (4.6),  the  other  with  (4.10).  Agreement  seems  quite 
satisfactory,  differences  appearing  at  most  in  the  fourth 
decimal  digit.  The  formula  (4.10)  needs  more  terms  to  be  a 


Table  6.  Comparison  of  Two  Formulae 

22  2  2  2  2 

z  =  (n-l)a£/a  a2/a 


800.0 

12.5017 

12.5017 

400.0 

8.5971 

8.5971 

200.0 

5.8367 

5.8367 

100.0 

3.8855 

3.8854 

50.0 

2.5069 

2.5067 

25.0 

1.5335 

1.5331 

12.5 

0.8475 

0.8467 

a2/a2  is  calculated  by  (4.5),  (4.6),  and 

a2/ a2  is  calculated  by  (4.10) . 

Si 
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good  approximation  for  smaller  values  of  z,  but  later  we 
will  see  that  the  case  with  z  less  than  25  is  not  important 
from  the  viewpoint  of  application. 

Note 

For  n  *  2,  the  formula  (4.10)  gives 

o^/a2  *  ~  a1/a'  -  f-  +  ^  a/o1  .  (4.11) 

In  Chapter  II,  the  exact  formula  of  the  optimal  ballistic 
dispersion  was  derived  for  n  =  2.  The  series  expansion  of 
a*/ a2  is 

cr^/a2  *  j  a^/a  -  -  +  oi/a ^  -  ...  (4.12) 

It  is  interesting  to  note  that  every  term  in  (4.11)  is  smaller 
than  the  corresponding  term  in  (4.12):  The  approximate 
value  given  by  (4.11)  for  n  »  2  is  always  smaller  than  the 
true  value. 

For  comparison's  sake,  we  will  also  consider  the  classi¬ 
cal  approximate  formula  by  Merritt  and  King  [Ref.  2]  who 
approximate  the  optimal  ballistic  dispersion  by  a 


where 


1  '  v'  T*T*7T TCP 


liuiaii  ii>|i 

.  *  .*  k  » 


ACCURACY  INVESTIGATION 


In  this  section  we  will  investigate  the  accuracy  of  both 

our  approximate  formula  and  the  classical  one. 

Table  7a  for  z  =  400,  Table  7B  for  z  =  100,  and  Table 

7C  for  z  =  25  illustrate  the  calculation  carried  out  for  the 

accuracy  study.  For  a  fixed  z,  the  approximate  value  for 

2 

the  optimal  ballistic  dispersion  divided  by  a  is  a  constant 

which  is  given  next  to  the  value  of  z  in  the  table.  On  the 

other  hand,  the  true  optimal  value  increases  with  n:  The 

first  column  is  n  and  in  the  second  column  is  given  the  corres- 

2  2 

ponding  true  optimal  a2/a  ,  the  calculation  of  which  is  based 

on  the  necessary  condition  of  optimality  (2.13),  and  requires 

fairly  long  calculation  time.  The  relative  error  of  our 

2  2  2  2 

approximate  formula  Act  =  (a  -a~)/a0  x  100%  is  given  in  the 

3rd  column.  Its  value  is  found  to  always  be  negative.  The 

salvo  kill  probability  associated  with  a2  *  0  is  denoted  as 

PQ  and  is  in  the  fourth  column.  The  fifth  column  is  the  upper 

bound  P  of  the  maximum  salvo  kill  probability  given  by  (4.6) 

and  (4.7)  .  In  the  sixth  column,  we  have  the  true  maximum 

salvo  kill  probability  designated  sintply  as  P  which  is  calcu- 

2 

lated  by  (2.21),  (2.22)  with  the  optimum  a2  given  xn  the 

2nd  column.  The  entries  in  the  7th  column  are  the  salvo  kill 

probabilities  associated  with  the  approximately  optimal  dis- 
2 

persion  a  .  The  notation  P  is  used  for  it.  The  last  column 

3  a 

is  the  relative  loss  in  the  kill  probability  which  would  be 
caused  if  we  were  to  use  the  approximate  value  a  instead  of 

cl 
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Table  7  Accuracy  of  Approximation  (4.10) 


z  *  400,  <r*/oc2  *  8.60 


n 

AO".1  \ 

P« 

5 

9.34 

-7.9 

0.0225 

10 

9.46 

-9.1 

0.0634 

15 

9.57 

-10.2 

0.1088 

20 

9.67 

-11.1 

0.1556 

25 

9.76 

-11.9 

0.2024 

30 

9.84 

-12.6 

0.2484 

35 

9.92 

-13.3 

0.2931 

40 

9.98 

-13.9 

0.3362 

nVoc* 

4.31 

4.41 

4.48 

4.55 

4.61 

4.66 

4.70 

4.74 


0.0414 

0.1484 

0.2777 

0.4063 

0.5229 

0.6230 

0.7061 

0.7733 


0.0410 

0.1456 

0.2709 

0.3952 

0.5081 

0.6056 

0.6873 

0.7542 


*  100  ,  Oi’/oc1  -3.89 


AO*.**/, 

-9.9 

-11.8 

-13.4 

-14.6 

-15.7 

-16.6 

-17.4 

-18.1 


0.0600 

0.1703 

0.2828 

0.3877 

0.4817 

0.5641 

0.6352 

0.6959 


0.0905 

0.2947 

0.4952 

0.6558 

0.7727 

0.8534 

0.9071 

0.9419 


0.0891 

0.2872 

0.4809 

0.6371 

0.7529 

0.8348 

0.8910 

0.9288 


*  25,  <rf(<xx  ■  1.53 


A  <*"„*  % 

Po 

1.77 

-13.6 

0.0573 

1.82 

-15.9 

0.2137 

1.86 

-17.8 

0.3706 

1.90 

-19.3 

0.5076 

1.93 

-20.5 

0.6206 

1.95 

-21.5 

0.7107 

1.98 

-22.4 

0.7812 

1.99 

-23.1 

0.8356 

0.0660 
0.2733 
0.4822 
0 . 6497 
0.7710 
0.8539 
0.9085 
0.9435 


0.0652 

0.2663 

0.4674 

0.6299 

0.7497 

0.8339 

0.8912 

0.9295 


0.0410 

0.1454 

0.2704 

0.3943 

0.5067 

0.6038 

0.6851 

0.7518 


0.0890 

0.2866 

0.4794 

0.6350 

0.7503 

0.8321 

0.8885 

0.9266 


0.0651 

0.2657 

0.4658 

0.6274 

0.7468 

0.8309 

0.8884 

0.9271 


2 

the  true  optimal  a  ,  AP  =  (P  -  p  )/P  *100%.  it  is  noted  that 

2  a 

in  the  case  with  a  large  z  value  the  salvo  kill  probability 

is  much  improved  by  employing  the  optimal  ballistic  disper¬ 
sion.  For  instance,  when  z  =  400,  the  ratio  P/Pg  is  around 

2.5.  But  the  improvement  is  not  so  much  when  z  gets  small. 

In  the  case  with  z  =  25,  the  ratio  is  at  most  1.25. 

As  to  our  approximation,  we  observe  the  following. 

2 

(1)  The  approximate  a  is  always  smaller  than  the  true 
value,  at  least  for  the  values  of  z  investigated. 

(2)  The  discrepancy  between  the  approximate  and  the 
true  values  is  larger  for  smaller  z. 

2 

(3)  For  a  fixed  z,  the  optimal  ballistic  dispersion  a2 
increases  as  n  is  increased.  Thus  the  discrepancy 
grows  when  n  increases  for  a  given  z. 

(4)  However,  we  may  say  that  the  discrepancy  is  not  so 
serious  because  the  kill  probability  associated  with 
the  approximate  optimal  ballistic  dispersion  is  not 
so  different  from  the  maximum  value.  The  relative 
loss  given  in  the  last  column  is  at  most  0.4%  which 
occurs  when  z  is  25,  a  less  important  case. 

Figure  5  is  presented  here  to  illustrate  the  properties 

2 

1  and  2.  The  ordinate  is  the  relative  error  Ao  ,  but  it 

cl 

2 

needs  several  words.  The  relative  error  A a  varies  with  n 

cl 

for  a  fixed  z  as  seen  in  Table  7,  To  get  a  representative 

2 

Aa  value  for  a  z,  an  n  is  picked  for  the  given  z  such  that 

61 

the  maximum  salvo  kill  probability  is  close  to  0.6,  an 


arbitrarily  chosen  figure.  Then  the  corresponding  Aa  is 
chosen  for  that  z  value,  for  instance  -12.63  for  z  =  400, 
and  similar  procedure  follows  for  other  z's.  The  bottom 
curve  in  Figure  5  shows  this  relationship.  Later  we  will 
discuss  other  curves  in  the  same  figure. 

From  the  figure  we  see  that  the  relative  error  is  always 
negative  and  its  absolute  value  is  decreasing  in  z.  The 
z-dependence  is,  however,  relatively  small,  and  we  learn 
its  bias  to  the  negative  side  is  a  characteristic  of  this 
approximate  formula. 

Similar  investigation  was  carried  out  on  the  classical 

approximate  formula  (4.13).  Table  8A  for  z  =  400,  Table  8B 

for  z  =  100,  and  Table  8C  for  z  =25  correspond  to  Table  7 
c  c 

the  only  difference  being  that  z  is  tabulated,  rather  than 

c 

z.  From  the  tables  we  observe  the  following. 

2  2 

(1)  The  approximate  a  /a  is  larger  than  the  true  value 

C 

except  for  the  cases  of  very  large  zc’s. 

(2)  The  discrepancy  between  the  approximate  and  the  true 

value  is  larger  for  smaller  This  tendency  is 

similar  to  that  of  our  approximation. 

.  .  2 

(3)  For  a  fixed  z ,  the  optimal  ballistic  dispersion  a2 
increases  as  n  is  increased.  Thus,  the  discrepancy 
gets  smaller  when  n  is  increased  for  a  given  zc  less 
than  800. 

(4)  However,  the  discrepancy  is  not  serious  because  the 
kill  probability  associated  with  the  approximate 
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Table  8  Accuracy  of  Approximation  (4.13) 


A  zc  =  400/  rel/ulm  10.00 


n 

5 

9.34 

7.1 

10 

9.46 

5.7 

15 

9.57 

4.5 

20 

9.67 

3.4 

25 

9.76 

2.5 

30 

9.84 

1.6 

35 

9.92 

0.9 

40 

9.98 

0.2 

n 

«\7a1 

4 

4.31 

15.9 

8 

4.41 

13.5 

12 

4.48 

11.5 

16 

4.55 

9.9 

20 

4.61 

8.5 

24 

4.66 

7.3 

28 

4.70 

6.3 

32 

4.74 

5.4 

C 


n 

Oa3/*1 

2 

1.77 

41.0 

4 

1.82 

37.1 

6 

1.86 

34.1 

8 

1.90 

31.7 

10 

1.93 

29.7 

12 

1.95 

28.0 

14 

1.98 

26.5 

16 

1.99 

25.3 

Po  P 

0.0225  0.0414 
0.0634  0.1484 
0.1088  0.2777 
0.1556  0.4063 
0.2024  0.5229 
0.2484  0.6230 
0.2931  0.7061 
0.3362  0.7733 


B  zc=  100,  tfVVo **•  5 


Po  P 

0.0600  0.0905 
0.1703  0.2947 
0.2828  0.4952 
0.3877  0.6558 
0.4817  0.7727 
0.5641  0.8534 
0.6352  0.9071 
0.6959  0.9419 


z,.3*  25,  c,1/**1  * 


Po  P 

0.0573  0.0660 

0.2137  0.2733 

0.3706  0.4822 

0.5076  0.6497 

0.6206  0.7710 

0.7107  0.8539 

0.7812  0.9085 

0.8356  0.9435 


P 

Pc 

*?% 

0.0410 

0.0410 

0.0 

0.1456 

0.1455 

0.0 

0.2709 

0.2709 

0.0 

0.3952 

0.3951 

0.0 

0.5081 

0.5080 

0.0 

0.6056 

0.6056 

0.0 

0.6873 

0.6873 

0.0 

0.7542 

0.7542 

0.0 

.00 


P 

Pc 

*P?o 

0.0891 

0.0889 

0.2 

0.2872 

0.2866 

0.2 

0.4809 

0.4800 

0.2 

0.6371 

0.6363 

0.1 

0.7529 

0.7523 

0.1 

0.8348 

0.8344 

0.1 

0.8910 

0.8907 

0.0 

0.9288 

0.9286 

0.0 

2.50 


P 

Pt 

4P2, 

0.0652 

0.0648 

0.5 

0.2663 

0.2638 

0.9 

0.4674 

0.4631 

0.9 

0.6299 

0.6249 

0.8 

0.7497 

0.7451 

0.6 

0.8339 

0.8300 

0.5 

0.8912 

0.8881 

0.3 

0.9295 

0.9272 

0.2 

optimal  dispersion  is  not  so  different  from  the  true 

maximum.  The  relative  loss  is  at  most  1%,  and  is 

2  2 

very  small  when  o^/a  is  large. 

2 

The  top  curve  in  Figure  5  is  the  relative  error  A a  of 

c 

the  classical  approximate  formula,  where  the  abscissa  should 
be  read  as  z  in  this  case.  In  spite  of  its  simple  form,  it 
is  indeed  an  excellent  approximation,  in  particular  for  large 
zc's.  The  zc“dependence  of  the  relative  error  A is,  how¬ 
ever,  more  sharp  than  ours.  The  middle  curve  in  Figure  5  will 
be  discussed  in  the  next  section. 

B.  IMPROVEMENT  OF  THE  APPROXIMATION 

In  the  preceding  section,  we  investigated  the  accuracy 
of  the  approximate  formula  given  by  (4.10),  and  it  seemed 
that  (4.10)  gave  a  lower  bound  to  the  true  optimal  ballistic 
dispersion.  The  reason  it  gives  a  possible  lower  bound,  and 
a  method  of  obtaining  improved  approximations  will  be  con¬ 
sidered  now. 

~  * 

The  Walsh  optimal  f2  and  others  in  (.3.16)  are  functions 
of  x, 

x  —  («J  +  u>2)  a2/2  ,  (4.15) 

for  the  Simple  Salvo  MOdel,  and  so  let  us  denote  them  as 
f2(x)  instead  of  f 2 (w^ ,w2) ,  etc.  With  this  notation,  (3.18) 
is 
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p*<X)/D(x)  . 


(4.16) 


^2 (x)  = 


2  2 

Equation  (3.27)  by  which  the  approximate  a  /a  is  determined 

3l 

reads 

2  .  2 

=  p0(x)/DQ(x)  ,  (4.17) 


f  2  ( x)  =  e 


-xaa/a 


and 


D(x)  =  DQ(x)  , 


(4.18) 


because  the  damage  function  is  assumed  to  be  Gaussian. 

The  function  p*(x)  is  the  double  Fourier  transform  of 
the  Walsh  optimal  p  function  given  in  Corollary  3.2  with 
°u  “  av  *  al '  and  Pq t*)  an  approximate  formula  for  it 
given  by  (3.26).  Figure  6  illustrates  what  the  curve  of  p*(x) 
is  like.  The  curve  decreases  as  x  is  increased  from  zero, 
changes  its  sign,  and  then  swings  back  to  positive  with  small 
amplitude,  and  so  forth.  On  the  other  hand,  the  approxi¬ 
mating  function  Pq Cx)  is  equal  to  p* (x)  at  x  =  0,  has  the 
same  tangent  there,  and  decreases  exponentially.  pQ (x)  is 
larger  than  p*(x)  for  a  wide  interval  of  x  as  is  seen  in 
Figure  6 . 


»  *  •  "  • 
»  *  .  • 


Wr< 
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m 
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pn(x)  >  p*(x) 


0  <  x  <  x. 
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Therefore,  from  (4.16)  and  (4.17),  we  have 

* 

e  >  f2(x)  ,  0  <  x  <  xQ 

2 

which  suggests  that  o  is  smaller  than  the  true  optimal 
value . 

One  of  the  simplest  ways  to  get  an  approximate  formula 

2 

would  be  to  define  a. 

o 

2  2 
-0X0, /a 

e  P0(x)/DQ(x)  (4.19) 

with  0  <  9  <  1,  or 

ab  *  aa/Q  '  0  <  6  <  1  ,  (4.20) 

and  to  determine  the  parameter  0  empirically  to  achieve 

.  2 
satisfactory  agreement  between  the  approximate  o^  and  the 

2  2 
true  a 2*  Note  that  the  definition  of  o^  corresponds  to 

approximating  the  p*(x)  function  by  p^(x)  with  an  exponent 

1/9  times  the  exponent  of  Pg (x) . 

For  example,  let  1/9  be  1.15.  Then, 

o£  =  1.15  a*  .  (4.21) 

In  Table  9A,  Table  9B,  and  Table  9C,  we  present  the  accuracy 
investigation  of  (4.21).  Accuracy  is  found  much  improved. 


Table  9  Accuracy  of  Approximation  (4.2/) 


A  z  =  400,  (T^/oi1  *  9.89 


n 

o~»  /«* 

A<rba% 

Po 

P 

P 

Pb 

5 

9.34 

5.9 

0.0225 

0.0414 

0.0410 

0.0410 

0.0 

10 

9.46 

4.5 

0.0634 

0.1484 

0.1456 

0.1455 

0.0 

15 

9.57 

3.3 

0.1088 

0.2777 

0.2709 

0.2709 

0.0 

20 

9.67 

2.3 

0.1556 

0.4063 

0.3952 

0.3952 

0.0 

25 

9.76 

1.3 

0.2024 

0.5229 

0.5081 

0.5081 

0.0 

30 

9.84 

0.5 

0.2484 

0.6230 

0.6056 

0.6056 

0.0 

35 

9.92 

-0.3 

0.2931 

0.7061 

0.6873 

0.6873 

0.0 

40 

9.98 

-1.0 

0.3362 

0.7733 

0.7542 

0.7542 

0.0 

B  z  »  100 

,  <r£  to*  * 

4.47 

n 

o-//*2 

4<rha% 

P. 

P 

P 

Pb 

4 

4.31 

3.6 

0.0600 

0.0905 

0.0891 

0.0691 

0.0 

8 

4.41 

1.4 

0.1703 

0.2947 

0.2872 

0.2872 

0.0 

12 

4.48 

-0.4 

0.2828 

0.4952 

0.4809 

0.4809 

0.0 

16 

4.55 

-1.8 

0.3877 

0.6558 

0.6371 

0.6371 

0.0 

20 

4.61 

-3.1 

0.4817 

0.7727 

0.7529 

0.7528 

0.0 

24 

4.66 

-4.1 

0.5641 

0.8534 

0.8348 

0.8347 

0.0 

28 

4.70 

-5.0 

0.6352 

0.9071 

0.8910 

0.8908 

0.0 

32 

4.74 

-5.8 

0.6959 

0.9419 

0.9288 

0.9286 

0.0 

C  z  *  25,  <rba/ 1.76 


n 

<*!*/«* 

4<rba7. 

P0 

P 

P 

Pb 

2 

-0.6 

0.0573 

0.0660 

0.0652 

0.0652 

0.0 

4 

1.82 

-3.3 

0.2131 

0.2733 

0.2663 

0.2663 

0.0 

6 

1.86 

-5.4 

0.3706 

0.4822 

0.4674 

0.4673 

0.0 

8 

1.90 

-7.1 

0.5076 

0.6497 

0.6299 

0.6296 

0.1 

10 

1.93 

-8.6 

0.6206 

0.7710 

0.7497 

0.7493 

0.1 

12 

1.95 

-9.7 

0.7107 

0.8539 

0.8339 

0.8333 

0.1 

14 

1.98 

-10.8 

0.7812 

0.9085 

0.8912 

0.8906 

0.1 

16 

1.99 

-11.6 

0.8356 

0.9435 

0.9295 

0.9239 

0.1 
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The  maximum  absolute  relative  error  A reduces  to  12%,  and 
the  relative  loss  in  the  salvo  kill  probability  is,  at  most 
only  0.07%.  The  middle  curve  in  Figure  5  is  the  relative 
error  A a? . 


V.  SALVO  MODEL  WITH  NON  GAUSSIAN  DAMAGE  FUNCTION 


In  this  chapter  we  deal  with  a  salvo  kill  model  with 
circular  normal  errors  and  circularly  symmetric  general 
damage  function.  Under  this  general  assumption,  calculation 
of  the  salvo  kill  probability  is  very  difficult,  and  time- 
consuming  computation  is  needed  for  it.  Here  we  propose  an 
approximate  method  for  calculating  the  salvo  kill  proba¬ 
bility.  By  this  method  computation  is  very  easy  and  yet 
fairly  high  accuracy  is  expected. 

Assumptions  in  this  chapter  are  as  follows:  A  salvo  of 
n  weapons  is  fired.  The  impact  point  error  is  composed  of 
two  parts,  the  random  bias  common  to  all  n  weapons  and  the 
round-to-round  error.  We  adopt  a  Cartesian  coordinate  sys¬ 
tem  (x,y)  the  origin  of  which  coincides  with  the  center  of 
impact  of  the  n  weapons.  The  random  position  of  the  target 
is  denoted  as  (U,V)  with  respect  to  this  coordinate  system, 

and  (U,V)  is  assumed  to  be  circular  normal  centered  at  C0,0> 
2 

with  variance  a^: 


f1(u,v)  =  — e  1  .  (5.1) 

2iTa^ 

The  impact  points  of  the  ith  weapon  are  independent 

and  identically  distributed  circular  normal  random  variables 


2 

centered  at  (0,0)  with  variance  a2: 


Ll3 


f2(x,y)  = 


2  2 
x  +y 

i 

- sc  e 


(5.2) 


The  conditional  kill  probability  of  the  target  by  a  single 
weapon  given  that  the  target  is  at  (u,v)  and  the  weapon  impact 
point  is  (x,y)  is  called  as  the  damage  function  D.  In  this 
chapter,  D  is  assumed  to  be  a  function  of  the  miss  distance 
r  =  >/(u-x)  2  +  (v-y) 2  only. 


D  =  D(r)  ,  r2  =  (u-x) 2  +  (v-y) 2  .  (5.3) 


Further,  the  lethal  area  is  denoted  as  A,  and  it  is  assumed 
that  there  exist  finite  ''moments"  of  all  degree  as  defined 
in  (5.5)  , 


2t r  /  D(r)  r  dr  =  A  ,  (5.4) 

0 


E(r^)  =  /  r^  D(r)  r  dr  <  00  ,  j  =  1,2,...  (5.51 

A  0 


The  conditional  kill  probability  by  a  single  weapon 
given  (.U,Vl  being  (u,v)  is  given  by  (.1.11,  repeated  here  as 


/  /  D (u-x, v-y) f, (x,y)  dx  dy  , 


(5 .6) 


—00  —CO 


p(u,v) 


and  the  salvo  kill  probability  (1.2)  is  renamed  as  (5.7). 


=  1-  /  /  (1-  p(u,v)  } 


n 


f  (u,v) 


du  dv 


(5.7) 


In  the  Simple  Salvo  Model  in  Chapter  II  with  Gaussian 
assumption  on  the  damage  function,  we  have  a  simple  expres¬ 
sion  for  p(u,v) : 


p  (u,  v)  = 


2 

a 

2';~2" 
a  +a2 


2  2 
u  +v 

2  2 
2  (ct  +a2) 


(2.6) 


In  the  general  model  here,  p(u,v)  given  by  (5.6)  is  a  com- 

5  2 

plicated  function  in  vu  +v  ,  and  computation  of  the  kill 
probability  is  difficult.  To  overcome  this  difficulty, 
we  propose  an  approximate  formula  for  p(u,v),  and  then 
present  a  recursive  method  for  computing  the  salvo  kill 
probability  with  this  approximate  p(u,v)  function. 

A.  APPROXIMATE  FORMULA  FOR  p(u,v) 

In  this  section,  an  approximate  formula  for  the  condi¬ 
tional  kill  probability  p(u,v)  is  presented.  As  a  corollary 
we  will  get  an  approximation  to  the  so-called  Circular 
Coverage  Function. 

Define  the  double  Fourier  transform  of  p(u,v)  as 
~  00  00  i(w,u+w2v) 

p  (u.fiiiJ  -  \  /  e  p(u,v)  du  dv  ,  (.5.8) 


with  similar  notations  being  used  for  D(x,y)  and  f(x,y). 
Then,  from  (5.6),  we  have 


p(a)1,(o2)  =  D  (oj^ /Ui2)  f2  *  (5.9) 

The  Fourier  transform  of  the  circular  normal  density  (5.2) 
is 


f 2  (a)i  /U)2} 


e 


2 

1 


) 


(5.10) 


In  the  Fourier  transform  of  the  damage  function,  let  the 
exponential  function  be  expanded  and  integrated  term  by  term, 


D  (oj^,(d2) 


»  oo  i(w1x+w,y) 

/  /  e  1  D(x,y)  dxdy 


00  00  oo  ,  • 

-If  /  Jf  iD  («1x  +u2y)  3D(x,y)  dxdy  . 

j  =  0  —co  —co  3  “ 


By  symmetry,  the  odd  terms  are  zero,  so 

®  k  .  oo  °°  . 

D(Wl,«2)  =  A  X  -gyj  J  /  /  C«1x+0)2y)  D(x,y)  dxdy 


-  A 


00 

X 


(-1) 

jm 


r  a 


ii 


i  2tt 

{/  (u. 


cos  0+w-  sin 


9) 2kd0}r2kD(r)r  dr 


Therefore, 


D  (u1  ,ui0)/h  =  1  +  T 
1  k=l 


(-l)k 

T2TcTT 


( 23c— 1 )  (2k-3)  •  •  -1 
2k(2k-2)  •  •  -2 


E(r 


2k\  ,  2.,  2. 


(5 


since 


2rr 


(2k-l) (2k-3) 
2kl2k-2 )  •  2“ 


k  ^ 


1  2k 

jjp  J  cos  0+0)2  sin  6)  <30  = 


k  - 


The  first  few  terms  in  (5.11)  are 


D(oj1,u2)/& 


1  -  J-  E(r2)  (u)2  +u2)  +  ^-E(r4)  (a)2  +  w2)2 


Now  let  us  define 


Da(u1/w2)  *  A{  1 


2  2 

h  2  2  "C (u. +W-) 

!(“i +“l),e  1  2 


(5 


and  determine  two  constants  b  and  c  such  that  (5.11)  and 
(5.12)  coincide  up  to  the  fourth  power  of  u^'s.  Then  we 
get 


/{E(r2}  }2  -E(r4)/2 


b 


1 

T 


i 


(5 


c  =  -b/2  +  E(r  )/4  . 


If  (5.10)  and  (5.13)  are  substituted  in  (5. 


pa(“l'w2) 


=  A{  1  - 


j(w2  +  u>2) }  e 


,  2\ 

■(ui^+o^) 


where 


S2  =  a2  -  b  +  E(r2)/2 

with  b  given  by  (5.13)  . 

Next,  the  inverse  transform  of  p  (u.  ,uj0) 

d  1  a 


Pa(u,v) 


1 

(  2tt  ) 


00 


7  / 

*-co 


-ifuw.+vw-)  - 
e  Pa(“! 


The  inverse  transform  of  the  first  term 


-S2(w2+w2)/2 
Ae  A  * 


is  easily  calculated  as 


u2+v 


2S‘ 


2ttS 


»)  ,  the  result  is 

,  (5.14) 

(5.15) 

is  calculated: 

(5.16) 

in  (5.14) 


The  second  term  in  (5.14)  is 


Ab  2  2  "  2  (wl+w2) 

— (co  i+to  2 )  e 


-b  ,  . 


d  ( S  ) 


and  its  inverse  transform  is  given  by 


2 ,  2 
u  +v 


^  d  r  A  _  2S , 

-b  - *— i - a-  e  } 

d(S^)  2ttS^ 


Therefore  we  have 


2  2 
u+v 


Pa (u,v) 

d 


2  2 

b  , ,  u  + v  \  i  „ 

(  1  "  )  J  ® 

S  2S^ 


(5.17) 


It  is  noted  that  (5.17)  is  exact  when  the  damage  function 

is  Gaussian  with  parameter  a;  (5.17)  coincides  with  the  exact 

2 

formula  (2.6)  because  in  this  case  A  =  2ira  ,  b  =  0,  and 

_2  2  ^  2 
S  —  o  2  +  cl  • 

The  cookie  cutter  damage  function  differs  substantially 
from  the  Gaussian,  so  consideration  of  the  former  provides 
a  good  test  of  the  robustness  of  (5.17)  .  The  cookie  cutter 
damage  function  is 


if  0  1  r  =  a 


D(r) 


(5.18) 


otherwise 


•■  M'.  1  V-  W  O 


The  probability 


00  oo 

p(u,v)  =  /  /  D (u-x, v-y)  f 2  (x,y)  dx  dy 

— OO  *00 


/  / 

x2+y2^a2 


(x-u) 2+ (y-v) 2 

2ol  ' 


dx  dy 


is  called  the  Circular  Coverage  Function.  In  the  literature, 

f~" 2  ? 

the  notation  P(a/a2,r/c2)  is  used  for  it,  where  r  =  vu  +  v  . 

A  standard  method  of  calculating  the  Circular  Coverage  Func¬ 
tion  is  the  algorithm  given  by  Brennan  and  Reed  [Ref.  10] 
based  on  an  infinite  series  expansion,  but  if  the  accuracy 
required  is  not  too  high,  it  is  convenient  to  use  closed  form 
approximations.  Standard  handbooks  [Refs.  11,12]  recommend 
the  next  formula  for  a/a2  <  1. 


PA(a/c2,r/c2) 


(5.19) 


with 


Z2  =  a2  +  a2/4  . 


(5.20) 


According  to  (5.17),  the  Circular  Coverage  Function 

(x-r) 2+y2 

1  2a2 

P(a/a2,r/a2)  =  J  J  .  *  dx  dy 


(5.21) 


2  x2+y2<a2 
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can  be  approximated  by 


Pa(a/a2,r/o2) 
where 

b  =  a2/4 /3  ,  s2  =  o2 +a2(l  -  l//3)/4  ,  (5.23) 

since  A  »  ira2,  E(r2)  =  a2/2,  and  E(r4)  =  a4/3  for  the 
cookie  cutter  damage  function. 

In  Table  10,  the  Circular  Coverage  Functio  P(R,d)  is 
tabulated  for  R  =  0.1  (0.1)  2.0  and  d  =  0.0  (0.5)  2.5.  In 
Table  11,  the  error  AA(R,d)  *  PA(R,dl  -  P(R,d)  is  shown. 

For  small  R,  the  error  is  very  small;  for  instance  the  error 
is  less  than  10  when  R  1  1.0.  But  when  R  is  increased 
beyond  1,  the  error  grows  rapidly. 

In  Table  12,  the  error  A  CR,d)  =  p  (R,d)  -  P(R,dl  of 

cL  a 

our  approximate  formula  is  given.  The  accuracy  is  found  much 

better  than  the  above-mentioned  formula.  The  absolute  error 

-4  <•  -2 

for  d  -  0  is  less  than  5  x 10  when  R  =  1.0,  less  than  10 

even  when  R  =  1.6  and  usually  shows  much  higher  accuracy. 

B.  RECURRENCE  FORMULA  FOR  THE  SALVO  KILL  PROBABILITY 

If  the  conditional  kill  probability  p(u,v)  by  a  single 
weapon  is  given  in  the  form  (5.17),  we  can  easily  compute  the 


rU  - 


2S4 


b  n  r 
—*•(1 - - 

S2  2S‘ 


■)  }e 


r 

2S‘ 


(5.22) 


Table  11  Error  Aa  (R#d) 


d 


R 

0.00 

0.50 

1.00 

0.1 

-0.0000 

-0.0000 

-0.0000 

0.2 

0.0000' 

0.0000 

0.0000 

0.3 

0.0000 

0.0000 

0.0000 

0.4 

0.0000 

0.0000 

0.0000 

0.5 

0.0001 

0.0001 

0.0000 

0.6 

0.0004 

0.0003 

0.0000 

0.7 

0.0010 

0.0007 

0.0001 

0.8 

0.0020 

0.0014 

0.0003 

0.9 

0.0038 

0.0027 

0.0005 

1.0 

0.0065 

0 . 0046 

0.0010 

1 . 1 

0.0106 

0.0076 

0.0018 

1 . 2 

0.0162 

0.0118 

0.0030 

1.3 

0.0236 

0.0174 

0.0048 

1 . 4 

0.0330 

0.0246 

0.0074 

1.5 

0.0447 

0.0337 

0.0109 

1.6 

0.0585 

0.0447 

0.0155 

1.7 

0.0746 

0.0578 

0.0213 

1.8 

0.0929 

0.0728 

0.0286 

1.9 

0.1132 

0.0898 

0.0374 

2.0 

0.1353 

0.1086 

0.0478 

1.50 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0001 

-0.0002 

-0.0004 

-0.0007 

-0.0012 

-0.0018 

-0.0026 

-0.0036 

-0.0047 

-0.0058 

-0.0068 

-0.0076 

-0.0080 

-0.0077 

-0.0065 


2.00 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0001 
-0.0002 
-0.0003 
—0 . 0006 
-0.0011 

-0.0019 

-0.0030 

-0.0045 

-0.0065 

-0.0090 

-0.0121 

-0.0157 

-0.0197 

-0.0240 

-0.0286 


■  v 


,v. 


to 


2.50 

•0 . 0000 
-0 . 0000 
■0.0000 
-0 . 0000 
-0.0000 

-0.0000 

-0.0000 

-0.0001 

-0.0002 

-0.0003 

-0.0007 

-0.0011 

-0.0019 

-0.0030 

-0.0046 

-0 . 0067 
-0.0095 
-0.0131 
-0.0174 
-0.0225 
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Table  12 


Error  Aa(R*d) 


d 


R 


0.00  0.50  1.00 


1.50  2.00  2.50 


0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 
1.9 
2.0 


-0.0000 

—0 . 0000 

—0 . 0000 

-0.0000 

—0 . 0000 

—0 . 0000 

-0.0000 

—0 . 0000 

—0 . 0000 

-0.0000 

—0 • 0000 

—0 . 0000 

-0.0000 

—0 • 0000 

0.0000 

-0.0000 

-0.0000 

0.0000 

-0.0000 

—0 • 0000 

0 . 0000 

—0 . 0001 

-0.0000 

0 . 0000 

-0.0001 

-0.0001 

0 . 0000 

—0 . 0003 

—0 . 0002 

0 . 0000 

—0 • 0006 

-0.0003 

0 . 0000 

-0.0011 

-0.0007 

0 . 0000 

-0.0019 

-0.0012 

0 . 0000 

-0.0031 

-0.0019 

0.0000 

-0.0049 

-0.0031 

—0 . 0000 

-0.0075 

-0.0048 

-0.0001 

-0.0110 

-0.0072 

-0.0003 

-0.0157 

-0.0104 

—0 . 0007 

-0.0218 

-0.0146 

-0.0012 

-0.0294 

-0.0199 

-0.0021 

— 0 . 0000 

—0 . 0000 

—0 . 0000 

— 0 . 0000 

-0 . 0000 

—0 . 0000 

— 0 . 0000 

-0.0000 

—0 • 0000 

0 . 0000 

—0 . 0000 

—0 . 0000 

0 . 0000 

0 . 0000 

—0 • 0000 

0 . 0000 

0 . 0000 

—0 . 0000 

0 . 0000 

0 . 0000 

—0 . 0000 

0 . 0000 

0 . 0000 

—0 . 0000 

0.0000 

0.0000 

—0 . 0000 

0.0001 

0 . 0000 

—0 . 0000 

0 . 0001 

0 . 0000 

—0 . 0000 

0 . 0003 

0 . 0001 

—0 . 0000 

0 . 0005 

0.0002 

-0.0001 

0.0008 

0.0003 

-0.0001 

0.0012 

0.0006 

-0.0001 

0.0019 

0 . 0009 

-0.0002 

0.0027 

0.0015 

-0.0002 

0.0039 

0.0023 

-0.0002 

0.0053 

0.0034 

—0 . 0001 

0.0071 

0.0049 

0 . 0001 

salvo  kill  probability  in  a  recursive  way.  We  first  observe 
that  the  salvo  kill  probability 


=  1“/  /  (1  -  p(u, v)  }nf1  (u,v)  du  dv 


—CO  —00 


where  f^(u,v)  is  the  circular  normal  density  (5.1)  and 


p(u,v)  -  .  (y  +  6  -£-y)  e  2S 

2S^ 


2  2  2 
,  r  =  u  +  v  , 


(5.24) 


can  be  expressed  as 


P  =  1  -  K  (n ,  0)  , 


(5.25) 


where 


oo  / 

K  (n,  0)  =  /  (1  -  (y  +  0t/p) e”t/p}ne~t  dt  , 

0 


15.26) 


2  2 

and  where  p  =  S  /a-^.  This  is  simply  a  matter  of  substitution 
and  the  introduction  of  polar  coordinates.  We  then  have 

Proposition  5.1 
The  function 


K(n,j)  =  /  (1  -  (y  +  St/p)e“t/p}ne“(j/p+1)  tdt 

0 


(5.271 
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r-.'  jr.y-r^T  -7. '7,'^*. V  -* ■;*. 


which  has  two  integer  arguments  and  three  parameters  y,  £ 
and  p  satisfies  the  following  recurrence  relation: 


K (n, j )  = 


n 


p - n_,,\n - Eli-  K(n-l,j+l)  , 


n+j+p  K(n_1'^  +  n+j+p (1_y)  “  n+j+p 

n  =  1,2,  ...  ;  j  =  0,1# •• • 


with 

K(0/ j)  = 

Proof 

From  (5.27) , 


.  ,  T  —  0/1/  . 

5  +  P 


(5.28) 


(5.29) 


i#r  1 


K(n, j)  -  /"{l- <y  +  6t/p)e-t/p}"'V(j/p+1)t  dt 
0 

0 

The  first  term  in  the  right-hand- side  is  K(n-1, j) .  The  second 
terms  is,  when  integrated  by  parts, 
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a£ai£ 


■* ZtLJkim,  .-M  *■-  -1  -* 


.V.  -  .'A 


-  /  U  -  <u  +Bt/p)e't/p>n‘1<u  +  St/p)e't/oe"<3/p+1)tdt 
0 

=  -  £<1-  (p  +  8t/c)e-t/p)n6-^/p+1)t!“ 

0 

-  £(j/p+l)  /  {1  -  (y  +  3t/p)e"t/p}ne"  <j/p+1)  t  dt 

0 

-  BrU-(u  +  0t/p)e-t/p}ne-t/pe-(=j/p+1)t  dt 

0 

=  £(l-u)n  -  i±£K(n,j)  -  BK(n-l, j+1)  . 

Therefore 

K(n,j)  =  K (n-1 ,  j )  +  £(l-y)  n  -  ^±2.  K(n,j)  -  f3K(n-l,j+l) 

from  which  (5.28)  is  immediate.  From  the  definition  of 
K(0,j),  it  is  clear  that  (5.29)  holds. 

Note  1 

We  use  the  notations  y  and  p  here  in  the  same  way  as 
in  Chapter  II  without  any  subscripts  because  no  confusion 
is  expected.  Their  definitions  are  different,  but  they  are 
just  the  counterparts  of  y  and  p  defined  in  Chapter  II. 

Note  2 

It  is  worthwhile  to  note  that  the  recurrence  relation 
(5.28)  has  a  similar  form  to  (2.17).  The  p  function  (5.24) 


involves  an  extra  parameter  3  which  was  missing  in  the  Simple 
Salvo  Model  in  Chapter  II,  and  therefore,  it  is  necessary 
to  introduce  another  index  j . 


VI.  COOKIE  CUTTER  SALVO  MODEL 


Consider  that  n  weapons  are  fired  at  a  target  in  a 
salvo.  The  impact  point  error  is  composed  of  two  parts, 
the  random  bais  and  the  dispersion.  We  adopt  a  Cartesian 
coordinate  system  (x,y)  the  origin  of  which  coincides  with 
the  center  of  impact  of  the  n  weapons.  The  random  position 
of  the  target  is  denoted  as  (U,V)  with  respect  to  this 

coordinate  system,  and  (U,V)  is  assumed  to  be  circular  normal 

2 

centered  at  (0,0)  with  variance  a^. 


fx (u,v) 


u2+v 


2 


(6.1) 


The  impact  points  of  the  ith  weapon  (X^Y^  are  indepen¬ 
dent  and  identically  distributed  circular  normal  random 

2 

variables  centered  at  (0,0)  with  variance 


f2 (**y) 


(6.2) 


The  kill  probability  by  a  single  weapon  conditional  on 
(U, V)  *  (u,v)  and  (X,Y)  *  (x,y)  is  a  function  of  u-x  and 
v-y,  D(u-x,v-y) ,  and  in  this  chapter  it  is  assumed  that 
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D  (u-x, v-y) 


(6.3) 


1  if  (u-x)2  +  (v-y)2  1  a2 
0  otherwise  , 

with  lethal  area 

_  2 
A  =  it  a 

In  the  following,  we  call  this  model  the  Cookie  Cutter 
Salvo  Model. 

The  conditional  probability  that  the  target  is  destroyed 
by  a  single  weapon  given  the  random  position  of  the  target 
being  (u,v)  is 

00  00 

p(u,v)  =  /  /  D (u-x, v-y) f 2 (x,y)  dx  dy  , 


with  (6.1)  and  (6.4)  is  readily  rewritten  as 


=  1-/{1  -  P (a/o 


2' 


^ne_t 


dt 


(6.5) 


We  carried  out  a  number  of  computations  of  (6.5)  with 
various  sets  of  parameters  c^,  a a,  and  n  by  the  Simpson 
method.  A  robust  algorithm  for  computing  the  Circular 
Coverage  Function  P(R,d)  for  wide  ranges  of  R  and  d  was  needed 
For  this  purpose,  we  used  the  Brennan- Reed  formula  in  a  modi¬ 
fied  form. 

Brennan- Reed  Formula: 


P(R,d) 


n=0 


9nkn 


9 


kn/kn-l  =  d2/ (2n)  ,  n  =  1,2,  ... 


-  (r2/2) n 


e'^/ru 


n  ■  1,2, 


Modified  Formula: 


(6.7) 


(6.8) 


The  modified  formula  is  especially  of  use  in  preventing 
underflow  possible  in  cases  with  large  R  and  d. 

In  their  review  paper,  Eckler  and  Burr  presented  a  short 
table  of  the  salvo  kill  probability  of  the  Cookie  Cutter 
Salvo  Model.  But  the  accuracy  is,  as  they  admitted,  sup¬ 
posedly  low  because  it  is  obtained  by  Monte  Carlo  simula¬ 
tion.  In  Table  13,  we  give  the  correct  values.  The  figures 
in  parentheses  are  from  the  table  by  Eckler  and  Burr. 

In  the  following  sections,  comparison  of  the  Cookie  Cutter 
Salvo  Model  and  the  Simple  Salvo  Model  is  dealt  with  and 
then  the  optimal  ballistic  dispersion  is  investigated,  where 
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Table  13  Salvo  Kill  Probability  :  Kisi/( Eckler-Burr ) 


A  n  38  2 


B  n  *  5 


C  n  *  20 


m 


/  O', 

0.25 

0.50 

0.75 

0.6 

0.221 

0.228 

0.199 

(0.224) 

(0.235) 

(0.201) 

1.2 

0.573 

0.587 

0.557 

(0.574) 

(0.582) 

(0.563) 

1.8 

0.834 

0.836 

0.813 

(0.834) 

(0.832) 

(0.814) 

2.4 

0.954 

0.952 

0.939 

(0.949) 

(0.951) 

(0.942) 

3.0 

0.991 

0.990 

0.985 

(0.990) 

(0.988) 

(0.985) 

Oi/cr, 

a/<r, 

0.50 

0.75 

1 .00 

0.4 

0.231 

0.213 

0.177 

(0.225) 

(0.214) 

(0.187) 

0.8 

0.525 

0.551 

0.514 

(0.521) 

(0.551) 

(0.516) 

1.2 

0.736 

0.772 

0.764 

(0.728) 

(0.773) 

(0.763) 

1.6 

0.872 

0.894 

0.896 

(0.865) 

(0.892) 

(0.893) 

2.0 

0.946 

0.957 

0.958 

(0.945) 

(0.956) 

(0.957) 

0.50 

a\/<ri 

0.75 

1.00 

1.50 

0.2 

0.221 

0.210 

0.176 

0.115 

(0.209) 

(0.201) 

(0.176) 

(0.113) 

0.4 

0.458 

0.525 

0.504 

0.381 

(0.447) 

(0.531) 

(0.499) 

(0.386) 

0.6 

0.605 

S.713 

0.740 

0.648 

(0.595) 

(0.716) 

(0.740) 

(0.640) 

0.8 

0.715 

0.818 

0.861 

0.828 

(0.706) 

(0.817) 

(0.846) 

(0.828) 

1.0 

0.801 

0.884 

0.922 

0.923 

(0.792) 

(0.882) 

(0.924) 

(0.924) 

1.2 

0.865 

0.928 

0.956 

0.966 

(0.856) 

(0.926) 

(0.958) 

(0.970) 

1.4 

0.912 

0.956 

0.975 

0.985 

(0.900) 

(0.953) 

(0.975) 

(0.984) 
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the  approximate  formula  presented  in  Chapter  III  and  the 
approximate  method  of  computing  the  salvo  kill  probability 
developed  in  Chapter  V  are  found  to  be  useful. 

A.  COOKIE  CUTTER  VS.  SIMPLE  SALVO  MODEL 

The  Cookie  Cutter  Salvo  Model  differs  from  the  Simple 
Salvo  Model  of  Chapter  II  only  in  the  damage  function.  In 
Chapter  II,  we  assumed 


(u-x)  2+  (v-y)  2 

- r  r  — 

D(u-x, v-y)  =  e  a  , 

2 

with  lethal  area  2tra  .  For  comparison's  sake,  the  lethal 
area  of  the  two  salvo  models  are  now  set  equal 

a2  =  2a2  .  (6.9) 

The  two  models  share  the  common  parameters  a^,  a a,  and  n, 
only  being  different  with  respect  to  the  shape  of  the  damage 
function.  The  quantities  associated  with  the  Cookie  Cutter 
Salvo  Model  will  be  given  suffix  c  in  the  following  when 
necessary . 

It  is  interesting  to  compare  the  salvo  kill  probabilities 
of  the  two  models.  We  begin  with  three  extreme  cases  in  which 
the  salvo  kill  probability  has  closed  form  expressions. 

(1)  The  case  n  =  1. 

The  formula  (6.5)  with  (6.4)  is  easily  integrated  and 


give 


P  =  1  -  e 

c 


-a  /(a1  +  a2) 


(6.10 


The  corresponding  single  shot  kill  probability  under  the 
assumption  of  the  Gaussian  damage  function  is 


P  =  1  - 


2  2  2 
1  +  a7(oJ  +  op 


(6.11 


Proposition  6.1 

The  Cookie  Cutter  Salvo  Model  with  n  =  1  gives  a  higher 
single  shot  kill  probability  than  that  of  the  Simple  Salvo 
Model. 


P  >  P  . 
c 


(6.12 


Proof 

Recall  that  e~x  <  l/(l+x)  for  x  >  0.  Let  a2/(a^  +  op  35 
in  (6.10)  and  (6.11),  then  (6.12)  is  obvious.  □ 


Generally,  in  the  salvo  with  n  >  1,  this  inequality 
does  not  necessarily  hold  as  will  be  shown  in  the  following 


(2)  The  case  =  0 
Let  a,  •*  0  in  (6.5)  ,  then  we  have 


The  corresponding  formula  in  the  Simple  Salvo  Model  is  given 
by  (2. 23),  namely 

P  =  1  -  1/ (1  +  a2/^) n  .  (6.14) 

(3)  The  case  =  0 

Let  a  j  0  in  (6.4)  .  We  get 

2,  2 

-a  /a, 

P  =  1  -  e  .  (6.15) 

c 


The  corresponding  formula  in  the  Simple  Salvo  Model  is  given 
by  (2.24)  : 


=  1  - 


_ n (n-1)  •  •  •  1 

(n+X)  (n-l+X)  •  •  • 


(6.16) 


with 


.  2.  2 

X  =  a  / a ^  . 

Proposition  6.2 

Between  the  salvo  kill  probabilities  of  the  Cookie 
Cutter  Salvo  Model  and  the  Simple  Salvo  Model  with  the  same 
parameters,  the  following  relationships  hold: 

If  s  0,  then 

P  >  P  .  (6.17) 

c 
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If  02  =  0,  then 

Pc  <  P  iff  o^/a2  <  1/XQ  ,  (6.18) 

where  XQ  is  a  positive  solution  to  the  equation 


n 

X  -  l  ln(l  +  X/j)  =  0  .  (6.19) 

j=l 


Proof 

The  proof  of  (6.17)  is  just  the  same  as  that  of  Proposi¬ 
tion  6.1.  As  to  the  second  assertion,  let 

P  =  1  -  e-9(X) 

in  (6.16) .  Then 


n  n 

g(X)  =  l  ln(l  +  X/j)  ,  dg/dX  =  £  l/(j+X)  . 

j=l  j-1 


Since  g(0)  =0,  dg/dX |  >  1, 

X=0 


2  2 
d  g/dX  < 


0 ,  we  have 


g(X)  <  X  iff  X  £  XQ  . 


where  XQ  is  given  by  (6.19)  .  The  relation  (6.18)  is 


immediate . 


□ 


As  to  the  general  case,  we  must  compare  the  two  kill 

probabilities  numerically .  To  illustrate  the  general  idea. 

Table  14  is  given.  In  Table  14A,  P  for  n  =  4  is  tabulated, 

and  the  corresponding  Simple  Salvo  Model  kill  probability 

P  is  in  Table  14 B.  It  is  observed  that  P  's  in  the  column 

c 

2  2 

a^/a  -  0  are  larger  than  P's  as  stated  in  Proposition  6.2. 

2  2 

P  's  in  the  first  row  corresponding  to  a,/ a  =  0  is  smaller 

C  e 

2  2 

than  P's  except  for  the  case  a^/a  =  0,  since  1/Aq  =  0.216 
for  n  =  4 . 

As  a  general  tendency,  P  is  larger  than  P  when  a_  is 

C  e 

relatively  larger  than  o^.  Pc  is  smaller  than  P  when  a 2 
is  considerably  smaller  than  o^,  as  is  seen  in  the  upper 
right  corner  of  the  table. 

B.  OPTIMAL  BALLISTIC  DISPERSION 

We  carried  out  a  golden  search  calculation  to  get  the 
2  2 

optimal  ratio  c^/a  for  a  number  of  cases  with  different 

2  2 

combinations  of  n  and  a^/ct  .  The  results  are  tabulated  m 

2  2 

Table  15.  Table  15A  gives  the  optimal  ^2c^a  t*ie  Cookie 

Cutter  Salvo  Model,  and  Table  15B  gives  the  associated 

salvo  kill  probability,  namely  the  maximum  salvo  kill  proba- 

2  2 

bility  for  a  given  set  of  n  and  a^/a  .  It  is  noted  that  this 
calculation  is  very  time-consuming.  The  corresponding 
values  for  the  Simple  Salvo  Model  are  tabulated  in  Table  16. 

Before  starting  comparison  of  the  optimal  ballistic 
dispersion  of  the  Cookie  Cutter  Salvo  Model  with  that  of  the 


Table  15  Optimal  Ballistic  Dispersion  and 
Maximum  Probability  (  Cookie  Cutter  Salvo  Model  ) 


A  Optimal  Ratio  /  <*  *■ 


2 

4 

8 

16 

32 

64 

2 

0.4 

1.0 

1.7 

2.7 

4.1 

6.1 

4 

0.8 

1.6 

2.6 

3.9 

5.8 

8.6 

6 

1.2 

2.1 

3.2 

4.9 

7.1 

10.5 

8 

1.6 

2.5 

3.8 

5.6 

8.3 

12.0 

10 

1.8 

2.9 

4.4 

6.4 

9.3 

13.4 

12 

2.1 

3.3 

4.9 

7.0 

10.1 

14.7 

B  Maximum  Salvo  Kill  Probability 

p 

*  C  max 


n 


2 

4 

8 

16 

32 

64 

2 

0.476 

0.626 

0.782 

0.905 

0.973 

0.996 

4 

0.300 

0.435 

0.601 

0.769 

0.899 

0.971 

6 

0.224 

0.339 

0.492 

0.667 

0.826 

0.935 

8 

0.182 

0.281 

0.420 

0.590 

0.762 

0.897 

10 

0.154 

0.241 

0.368 

0.531 

0.708 

0.859 

12 

0.135 

0.213 

0.329 

0.484 

0.661 

0.823 
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Table  16  Optimal  Ballistic  Dispersion  and 


Maxmum  Probability  (  Simple 


A  Optimal  Ratio 


2 

4 

8 

0.5 

0.6 

1.3 

0.5 

1.1 

2.0 

0.5 

1.4 

2.6 

0.7 

1.8 

3.1 

0.9 

2.1 

3.6 

1.0 

2.3 

4.0 

B  Maximum  Salvo  Kill 

^  «n*  x 


2 

4 

8 

0.450 

0.622 

0.780 

0.293 

0.431 

0.599 

0.217 

0.334 

0.489 

0.173 

0.274 

0.416 

0.144 

0.233 

0.363 

0.124 

0.204 

0.323 

Salvo  Model  ) 


ff»  /d1  . 


16 

32 

64 

2.2 

3.6 

5.6 

3.4 

5.3 

8.1 

4.3 

6.6 

10.0 

5.0 

7.7 

11.5 

5.7 

8.7 

12.9 

6.3 

9.5 

14.1 

Probability 


16 

32 

64 

0.904 

0.973 

0.996 

0.768 

0.899 

0.971 

0.665 

0.826 

0.935 

0.588 

0.761 

0.896 

0.528 

0.707 

0.858 

0.481 

0.659- 

0.823 

Now  we  subtract  entities  in  Table  16A  from  those  in  Table 

2  2  2  2 

15A  and  get  the  difference  cr‘c/a  -  which  tabulated  in 

Table  17a.  The  values  are  positive  with  only  one  exception, 
and  the  predicted  value  0.5  is  observed  in  the  upper  right 
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Table  17  Comparison  of  Cookie  Cutter  Salvo  Model 

and  Simple  Salvo  Model 


A  Difference  of  Optimal  Ratio  <Txc /<xx  -  Gl'/oC*1 


n 


Oiid 

2 

4 

8 

16 

32 

64 

2 

-0.1 

0.4 

0.4 

0.5 

0.5 

0.5 

4 

0.3 

0.5 

0.6 

0.5 

0.5 

0.5 

6 

0.7 

0.7 

0.6 

0.6 

0.5 

0.5 

8 

0.9 

0.7 

0.7 

0.6 

0.6 

0.5 

10 

0.9 

0.8 

0.8 

0.7 

0.6 

0.5 

12 

1.1 

1.0 

0.9 

0.7 

0.6 

0.6 

B 

Difference  of 

Maximum  P 1 s 

Pcm«j* 

P 

m  «* 

n 

<**/«* 

2 

4 

8 

16 

32 

64 

2 

0.026 

0.004 

0.002 

0.001 

0.000 

0.000 

4 

0.007 

0.004 

0.002 

0.001 

0 .000 

0.000 

6 

0.007 

0.005 

0.003 

0.002 

0.000 

C.001 

8 

0.009 

0.007 

0.004 

0.002 

0.001 

0.001 

10 

0.010 

0.008 

0.005 

0.003 

0.001 

0.001 

12 

0.011 

0.009 

0.006 

0.003 

0.002 

C.001 
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corner  of  the  table.  In  the  below  left  corner  of  the  table, 

the  figures  are  not  close  to  the  value  0.5,  but  these  cases 

correspond  to  low  salvo  kill  probability. 

In  Table  17B,  the  difference  of  the  maximum  probabilities 

of  each  model  is  tabulated.  The  figures  are  all  positive, 

and  lead  to  the  speculation  that  P„ _  „  >  P _ _  in  all  cases 

c  max  max 

despite  the  fact  that  Pc  may  be  smaller  than  P  when  a 2  is 
not  chosen  optimally.  The  author  has  not  yet  succeeded  in 
obtaining  a  proof  of  this  speculation,  but  it  might  be  the 
counterpart  for  n  >  1  of  Proposition  6.1  for  n  *  1. 

In  Chapter  V,  we  developed  an  approximate  method  of 
computing  the  salvo  kill  probability  which  can  be  applied 
to  the  Cookie  Cutter  Salvo  Model.  The  salvo  kill  proba¬ 
bility  is  given  by 


P  =  1  -  K(n, 0}  , 


(5.25). 


where 


K(n,j)  -  /  { 1-  (y  + Bt/p)e-t/p}ne'(j/p+1)t  dt  (5.27) 

0 


satisfies  the  recurrence  relation. 


ferKU-l.j)  H-=rSrra-u)n-=^rK(n-l,j+l),  (5.28) 


K(n, j) 


with 


K(0,j)  - 


5+p  ' 


j  =  1,2, 


(5.29) 


The  parameters  are  given  by  (5.22),  (5.23),  and  (5.24). 


y  =  (a2/S2)  (1  -  a2/2/3  S2)  , 


8  =  a4/2/3  S4  , 


(6.22) 


S2  =  a2  +  a2  (1  -  l//J)/2 


Using  these  formulae  together  with  the  golden  section 
search  calculation,  we  obtained  Table  18.  Table  18A  gives 
the  approximate  optimal  ballistic  dispersion  of  the  Cookie 
Cutter  Salvo  Model,  and  the  associated  salvo  kill  probability 
is  in  Table  18B. 

2  2 

In  Table  19A,  the  difference  between  the  optimal  a2c/a 

2  2 

given  in  Table  15a  and  the  approximate  c^a^01  Tab^e  ^8A 
is  tabulated.  The  difference  decreases  towards  the  upper 
right  corner  of  the  table. 

In  Table  19 B,  the  salvo  kill  probability  of  the  Cookie 
Cutter  Salvo  Model  associated  with  the  approximately  optimal 
ballistic  dispersion  mentioned  above  is  given.  It  corres¬ 
ponds  to  the  maximum  P  given  in  Table  15B.  In  spite  of  the 

2  2 

discrepancies  observed  in  Table  19A,  the  approximate  w* 

2  2 

gives  almost  the  same  kill  probability  as  the  true  °2c'a  91veS- 
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Table  18  Optimal  Ballistic  Dispersion 
and  the  Associated  Salvo  Kill  Probability 
(  Cookie  Cutter  Salvo  Model,  Approximation  ) 


o-*/otx 

A 

2 

Optimal 

4 

.  Ratio 

rv 

8 

CzA*  /<XJ  • 

16 

32 

64 

2 

0.5 

1.0 

1.7 

2.7 

4.1 

6.1 

4 

0.7 

1.5 

2.5 

3.9 

5.8 

8.6 

6 

1.0 

1.9 

3.1 

4.8 

7.1 

10.5 

8 

1.2 

2.2 

3.6 

5.5 

8.2 

12.0 

10 

1.3 

2.5 

4.1 

6.2 

9.2 

13.4 

12 

1.5 

2.8 

4.5 

6.8 

10.0 

14.7 

B 

Maximum 

Salvo  Kill 

Probability 

*> 

n 

O-.'/ot* 

2 

4 

8 

16 

32 

64 

2 

0.475 

0.625 

0.782 

0.905 

0.973 

0.996 

4 

0.297 

0.432 

0.600 

0.768 

0.899 

0.971 

6 

0.218 

0.334 

0.490 

0.665 

0.826 

0.935 

8 

0.173 

0.274 

0.416 

0.588 

0.762 

0.396 

10 

0.144 

0.234 

0.363 

0.528 

0.707 

0.358 

12 

0.124 

0.204 

0.323 

0.481 

0.659 

0.323 
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Table  19  Accuracy  of  the  Approximation 
{  Cookie  Cutter  Salvo  Model  ) 


A  Error  in  the  Optimal  Ratio,  -  <r2A2  /o(4  . 


n 


2 

4 

8 

16 

32 

64 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.1 

0.1 

0.1 

0.0 

0.0 

0.0 

0.2 

0.2 

0.1 

0.1 

0.0 

0.0 

0.4 

0.3 

0.2 

0.1 

0.1 

0.0 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 

0.6 

0.5 

0.4 

0.2 

0.1 

0.0 

B  Salvo  Kill  Probability  Associated  with 
the  Approximately  Optimal  Ballistic  Dispersion. 

?  c  <r»Aa) 


n 


7otX 

2 

4 

8 

16 

32 

64 

2 

0.475 

0.626 

0.782 

0.905 

0.973 

0.996 

4 

0.300 

0.435 

0.601 

0.769 

0.899 

0.971 

6 

0.223 

0.339 

0.492 

0.667 

0.826 

0.935 

8 

0.181 

0.280 

0.420 

0.590 

0.762 

0.397 

0 

0.153 

0.241 

0.368 

0.531 

0.708 

0.359 

2 

0.134 

0.212 

0.328 

0.484 

0.661 

0.323 
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VII.  SALVO  MODELS  WITHOUT  CIRCULAR  SYMMETRY 


The  approximate  formula  for  the  optimal  ballistic  dis¬ 
persion  presented  in  Chapter  III  covers  the  salvo  models 
without  circular  symmetry.  In  this  chapter  the  formulae 
for  computing  the  salvo  kill  probability  of  these  models 
are  dealt  with. 

Suppose  that  a  salvo  of  n  weapons  is  fired  against  a 

target.  The  impact  point  error  is  composed  of  two  parts, 

one  being  common  to  all  the  weapons,  and  another  being 

round-to-round  dispersion.  The  random  target  position  (U,V) 

with  respect  to  the  coordinate  system  (x,y)  is  assumed  to 

be  elliptical  normal  centered  at  (0,0)  with  variances 
2  2 

<tu  and  av  respectively. 


The  kill  probability  D  by  a  single  weapon  conditional  on  the 

target  being  at  (u,v)  and  the  weapon  impact  point  being  (x,y) 

is  a  function  of  u-x  and  v-y  only,  and  is  given  by  the 

elliptical  Gaussian  function  with  parameters  a  and  a  . 

x  y 

Specifically, 


D(u-x, v-y) 


(u-x) ' 
x 


(v-y) 

2a 


(7.3) 


It  might  be  an  appropriate  model,  e.g.,  for  a  ground  target 
vs.  a  weapon  impacting  the  ground  obliquely,  because  in  this 
case  the  scattering  of  the  fragments  and  thus  the  conditional 
kill  probability  D  is  by  no  means  circular  symmetric  as 
assumed  in  the  previous  chapters.  For  brevity,  we  call 
this  model  the  Elliptical  Normal  Salvo  Model. 

The  kill  probability  of  a  single  weapon  conditional  on 
(U, V)  being  (u,v)  is  given  by 


CO  CO 

p(u,v)  =  /  /  D(u-x, v-y)  fj  (x,y)  dxdy  , 

—00  —00 


and  using  (7.2)  and  (7. 


we  get 


p(u,v) 


ctxcty 


u 

T 


2  «4+0J) 


(7.4) 
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The  salvo  kill  probability  is  then  given  by 


P 


1  - 


1 

2  it  a  a- 

U  V 


/ 


/  {l-p(u,v) }n 


—00 


v  dudv  . 


(7.5) 


From  this  equation,  Grubbs  derived  the  following  formula: 
Proposition  7.1  (Grubbs) 

The  kill  probability  of  the  Elliptical  Normal  Salvo 
Model  stated  above  is  given  by 


N 


P  = 


l  (-D 
j-1 


j-1 


pxp 


(d^T 


3+p^" 


(7.6) 


where 


v  “  axay/  v/<0tx  +  ax}  <ay  +ay)  ' 

p  *  (a 2+a2)/a2  ,  p  =  (a2+a2)/o2  . 

x  x  x  u  y  y  y  v 


(7.7) 


The  Grubbs  formula  is  not  suited  for  calculation  when  n 
is  large,  since  it  is  an  alternating  series.  To  overcome  this 
difficulty,  Breaux  and  Mohler  [Ref.  13]  gave  a  method  for  calcu¬ 
lating  the  kill  probability  based  on  an  expansion  of  (l-z)n 
in  Jacobi  polynomials  rather  than  as  a  binomial  series.  The 
series  is  found  to  converge  to  the  true  value  with  less  than 
n  terms,  which  is  very  attractive  for  the  calculation  of 
salvo  kill  probability  with  a  large  n. 


The  author  has  tried  to  find  alternative  ways  of  calcu¬ 
lating  the  salvo  kill  probability  which  are  effective  even 
for  the  cases  with  large  n's,  but  at  present,  the  obtained 
results  are  not  promising.  The  following  propositions  are 
given  not  as  efficient  algorithms  for  calculating  the  salvo 
kill  probability,  but  as  possible  hints  for  developing  more 
efficient  ways  of  calculation. 

For  the  Elliptical  Normal  Salvo  Model,  we  have  a 
proposition  similar  to  Proposition  2.5.  We  omit  the  proof, 
since  it  closely  parallels  that  of  Proposition  2.5. 


Proposition  7.2 

The  salvo  miss  probability  Q(n)  of  the  Elliptical  Normal 
Salvo  Model  characterized  by  (7.1) ,  (7.2) ,  and  (7.3)  with 
the  weapon  number  n. 


Q(n) 


/  (1  -  ye' 


2  «-2 
•s  -t 


n 


-P. 


e 


2  J 

s  -P  t 

Ky 


ds  dt 


(7.8) 


satisfies  the  relationship 


Q(n)  «  l  (?)  (l-y)n"k  yk  q(k)  , 


(7.9) 


k-0 


where 


q(k) 


,/o  0  00  00  2  ,  2  _  -p  s^”P 

/  /  u-e-»  )n.  *  y 

"  —00  -00 


y  ds  dt  .  (7.10) 
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In  the  case  of  Simple  Salvo  Model,  q(k)  had  a  simple 


recurrence  relationship,  but  we  have  not  succeeded  in 
finding  such  a  relationship  in  this  case. 

Corresponding  to  the  recurrence  formula  for  the  salvo 
kill  probability  presented  in  Chapter  II,  we  have  the  next 
propositions:  The  series  given  in  Proposition  7.3  is  an 
infinite  series,  but  for  a  small  6/p  it  will  converge  rapidly 

Proposition  7.3 

The  salvo  miss  probability  of  the  Elliptical  Normal 
Salvo  Model  with  n  weapons  is  given  by 


Q(n) 


A  -62/pi  l  ^2k)  S  (S/p)2k  Jtn,2k)  , 
k=0  2K(kir 


(7.11) 


where 


J(n,k) 


/  (1  -ye“t/p)n  e-t  tk  dt  , 


FT 


with 


2  .  2 

T" 

T 


2  2 
a..  +  a. 


'u 


2  2 
l,ax  +  0x 
T - j— 


a2  +  a2 


(7.12) 
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I*  N**.’  *  V*  •  •  * 


and 


y  ■ 


ct  a 

x.y 


y/  (a2 


K  +  «*>  «°y  +  °y> 


Proof 

In  (7.8),  let  s  =  r  cos  0 ,  t  =  r  sin  6 ,  and  (px+p^)/2 
(px-py)/2  =  fi.  Then 


W  0  "  m 


2  2  2 
p  s  +  p  t  =  r  (p  +  5  cos  28) 
x  y 


[A  A' 


and  therefore  Q(n)  is 


Q(n)  *  — /  (l-ue"r  )ne"pr2{/2lT  e"6r2  cos  20de}r  dr 

0  0 


<»  2  2 

2,/Tir  /  (1  -  ye"r  )  ne”pr  I  (6r^)  r  dr  , 

*  y  0  u 


where  IQ(z)  is  the  modified  Bessel  function  of  the  first 
kind  with  order  0,  and  has  an  expansion  formula 


v*> 


1  1  2 ,z, 2k 


I  %)  n|) 

k=0  Kl  * 


(7. 


Thus 


p » 
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Q (n)  =  /pxPy/P  I  a-ue“t/p)ne“tl0(6t/p)dt 


t  r~  r  (2k)!  /r/.v^K  1  r  ..  _-t/pv  n.-t^jc^ 

pxpv  He - 5‘(|5/P)  TT7TT  j  )  «  t  dt 

x  y  k=0  2*(k!)  *  UJc;  *  0 


2k  1 


-t/p. n  -t^2k. 


Proposition  7.4 


The  integral  J (n,k)  which  has  two  integer  arguments  and 


two  positive  parameters  u  and  p. 


J(n,k) 


•t/ p  n  -t  .  k 


)“  e  t  dt 


satisfies  the  following  recurrence  relation. 


J(n,k)  = 


Sp  +  H^p  J'n'k-1>  - 


n  =  1,2,  . . .  ;  k  =  1,2, 


J (n,0) 


-  n?p  J(n-1,0>  ♦  (1-U)n  ,  n  -  1,2, 


J(0 ,k)  *  1  ,  k  ■  1,2,  ... 


(7.14) 


(7.15) 


(7.16). 


(7.17) 


Proof 


J(n,k)  =  pr  /  (1  -ye’^)11"1  e_t  tk  dt 

K.  q 

1  r°°/i  -t/p .  n-1  -t/p-t  .k 
-  p  rr  i  (1  -ye  )  e  t  dt  . 

0 

The  first  term  in  the  right-hand-side  is  J(n-l,k).  If 
the  second  term  is  integrated  by  parts,  then 

00 

J(n,k)  =  J(n-l,k)  -  ^  £  (1 -ue“t/p)ne-t  tk| 

+  |  /q(  1  -  ye^^P)  ne-t  k  t1*"1  dt 

-  ITT  &  /”(1  *  “e‘t/P>  "•‘t  tk  at 

Therefore,  for  k  f  0,  we  have 
J (n,k)  *  J (n-l,k)  +  £  J(n,k-1)  -  £  J(n,k)  . 

From  this  relation  (7.15)  is  immediate.  If  k  =  0, 

J(n,0)  -  J(n-1,0)  +  £  (1  -p)n  -  £  J(n,0)  , 

and  therefore,  (7.16)  results.  The  relation  (7.17)  is  obvious. 

□ 
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It  is  noted  that  all  the  terms  in  the  expansion  (7.11) 
are  positive,  and  that  the  integral  J(n,0)  is  the  salvo  miss 
probability  of  the  Simple  Salvo  Model  dealt  with  in  Chapter 


II.  Therefore  we  have 
Corollary  7.5 

A  lower  bound  to  the  salvo  miss  probability  of  the 
Elliptic  Normal  Salvo  Model  characterized  by  (7.1),  (7.2) 
and  (7.3)  with  weapon  number  n  is  given  by 


QE(n) 


Q  (n)  , 


(7.18) 


where  Q(n)  is  the  miss  probability  of  the  Simple  Salvo  Model 
by  n  weapons  with 


THtJ 
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